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ABSTRACT 


This  thesis  demonstrates  through  three  comprehensive  examples, 
the  capabilities  of  the  TI  59  programmable  hand-held  calculator  as  an 
analytical  tool.  One  example  is  a  probability  application  while  the 
other  two  examples  entail  use  of  the  TI  59  in  statistical  inference  and 
data  analysis.  The  probability  example  involves  the  use  of  the  Monte 
Carlo  technique  to  simulate  stochastically  the  detection,  identification 
and  engagement  of  a  cruise  missile  by  an  Improved  Hawk  Air  Defense 
Battery.  The  second  example  illustrates  a  TI  59  program  which  is  de¬ 
signed  to  analyze  sample  data.  The  data  used  for  this  illustration  were 
gathered  by  the  authors  in  an  experiment  which  encompassed  the  testing 
of  thirty-six  male  subjects  to  determine  the  extent  to  which  their 
training  routines  influenced  their  strength,  endurance,  and  cardio¬ 
vascular  fitness.  The  third  example  involves  the  use  of  an  ANOVA  routine 
and  Scheffe's  multiple  contrasts  to  demonstrate  how  the  TI  59  may  be  used 
to  facilitate  statistical  inferences.  The  fitness  data  are  also  used 
for  this  purpose.  The  intent  throughout  the  thesis  is  to  exemplify  the 
capabilities  of  the  TI  59  as  a  viable ,  real  world  analytical  tool  rather 
than  emphasize  particular  results  of  the  simulation  or  the  experiment. 
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I .  INTRODUCTION 


The  intent  of  this  thesis  is  do  demonstrate  through  three  compre¬ 
hensive  examples,  the  tremendous  capabilities  of  the  TI  59  programmable 
hand-held  calculator.  One  of  the  examples  is  a  probability  application 
while  the  other  two  entail  use  of  the  TI  59  in  statistical  inference 
and  data  analysis. 

The  example  chosen  to  illustrate  an  application  to  probability 
theory  is  a  combat  simulation  model.  The  model  involves  use  of  the 
Monte  Carlo  technique  to  simulate  stochastically  detection,  identifi¬ 
cation,  and  engagement  of  a  cruise  missile  by  an  Improved  Hawk  Air 
Defense  Battery.  Chapter  III  discusses  the  combat  model  and  the  TI  59 
simulation  in  detail.  Chapter  II  addresses  briefly  the  Monte  Carlo 
technique,  combat  modeling  in  general,  and  the  TI  59  random  number 
generator  in  order  that  the  reader  may  better  understand  the  combat 
model  discussed  in  Chapter  III. 

Two  examples  are  provided  for  statistical  applications  of  the 
calculator.  These  involve  the  analysis  of  data  gathered  by  the  authors 
in  a  physical  fitness  experiment.  The  experiment,  discussed  in  detail 
in  Chapter  IV,  involved  the  testing  of  thirty-six  male  subjects  who  were 
divided  into  six  categories  based  on  their  training  routine,  to  deter¬ 
mine  whether  the  subjects'  training  program  did  in  fact  influence  their 
physical  fitness.  Chapter  IV  describes  the  fitness  experiment  in  terms 
of  its  scope,  design  and  scoring  methodology.  Chapter  IV  also  discusses 
a  TI  59  program  which  computes  measures  of  central  tendency  and  spread 
and  then  illustrates  the  use  of  the  program  with  the  fitness  data. 


Chapter  V  describes  a  TI  59  program  for  analysis  of  variance  and  then 
demonstrates  how  the  program  may  be  used  with  fitness  data  to 
make  statistical  inferences. 


Throughout  the  thesis,  it  is  assumed  that  the  reader  is  generally 
familiar  with  programming  techniques  for  the  Tl  59  hand-held  calculator. 
Subroutines,  labels,  flags,  data  registers,  and  program  steps  are 
discussed  in  each  of  the  ensuing  chapters  where  the  intent  is  to 
illustrate  how  the  features  of  the  calculator  may  be  exploited  to 
facilitate  statistical  analysis  or  simulation.  Reference  12  discusses 
programming  techniques  for  the  TI  59 .  The  TI  59  has  one  particular 
feature  which  makes  it  much  more  than  a  calculator.  Specifically  the 
capacity  to  use  subroutines  provides  a  analytical  tool  more  like  a 
minicomputer  than  a  calculator.  The  three  programs  discussed  in  the 
succeeding  chapters  use  subroutines  extensively  to  illustrate  this 
powerful  capability. 

Finally,  while  a  few  of  the  referenced  tables  and  charts  of  this 
thesis  are  positioned  close  to  comments  discussing  their  purpose, 
most  are  to  be  found  at  the  conclusion  of  each  chapter  or  in  the 
appendices. 


Q 


II.  PROBABILITY 


The  intention  of  the  authors  was  to  begin  this  chapter  discussion 
with  a  definition  of  probability  theory,  that  branch  of  mathematics 
generally  believed  to  have  been  founded  by  a  Swiss  mathematician  named 
Jacques  Bernoulli.  However,  research  has  revealed  that  there  is  some 
discussion  as  to  the  true  meaning  of  probability  theory  and  that  among 
mathematicians  there  appear  to  be  those  who  view  probability  as  a  state 
of  the  universe  while  others  consider  probability  a  state  of  belief. 

To  compound  this  situation  furthermore ,  there  appear  to  be  differing 
definitions  of  probability  within  each  group.  Indeed,  all  attempts  to 
define  probability  directly  have  failed  to  meet  with  success.  Instead, 
probability  has  been  axiomitized,  much  like  geometry,  so  that  a  set  of 
consistency  rules  or  axioms  established  by  A.  Kolmogorov  are  now 
generally  accepted  by  the  scientific  community.  These  axioms  allow  a 
great  deal  of  freedom  in  the  assignment  of  probabilities  for  any  par¬ 
ticular  model  a.xd  at  the  same  time  force  any  such  assignment  to  be 
consistent  with  any  other.  Moreover,  the  theorems  of  that  theory  then 
become  universally  true  statements  for  any  such  assignment.  In  this 
system,  events  are  defined  as  sets  in  a  specified  sample  space.  With 
those  guidelines  as  a  background,  probability  theory  can  be  used  to  make 
intelligent  predictions  and  decisions  if  we  know  what  events  are  possible 
and  how  probable  are  the  various  events.  After  a  little  thought  it  be¬ 
comes  immediately  apparent  that  the  immense  power  of  such  a  tool  as 
probability  theory  is  limited  in  use  only  by  one’s  imagination  and 
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ingenuity.  This  research  is  u.  effort  to  use  probability  theory  in  the 
construction  of  a  probabilistic  combat  simulation  on  the  Texas  Instruments 
programmable  59  calculator  (TI  59).  Because  the  simulation  developed 
includes  a  number  of  the  many  chance  elements  involved  in  most  combat 
situations,  a  discussion  of  the  Monte  Carlo  technique  and  random  number 
generation  on  the  TI  59  follows.  A  brief  disucssion  of  combat  model 
simulations  concludes  this  chapter. 


A.  MONTE  CARLO  TECHNIQUE 


Systems  that  exhibit  stochastic  elements  in  their  behavior  can  be 
simulated  with  the  aid  of  the  technique  called  Monte  Carlo  (named  after 
the  famous  gambling  resort  town  of  Monaco).  This  technique  involves 
sampling  from  those  known  probability  distributions  that  represent  each 
of  the  actual  chance  processes  included  in  the  system  under  study  [Bef.  3 
By  completing  a  system  simulation  run  many  times  while  keeping  the  non¬ 


stochastic  inputs  constant  but  allowing  the  chance  elements  to  fluctuate 
according  to  their  known  probability  distributions,  a  statistical 


average  for  run  results  can  be  determined. 


Turban  and  Meredith 


Ref.  8 


have  listed  the  steps  necessary  in 


building  a  Monte  Carlo  simulation  as  follows : 


"1.  Describe  the  system  and  obtain  the  probability  distributions 
of  the  relevant  probabilistic  elements  of  the  system. 

2.  Define  the  appropriate  measure (s)  of  performance. 

3.  Construct  cumulative  probability  distributions  for  each  of 
the  stochastic  elements. 

4.  Assign  representative  numbers  in  correspondence  with  the 
cumulative  probability  distributions. 

5.  Generate  a  random  number  for  each  of  the  independent  stochastic 
elements  and  .  .  .  (determine)  the  measure  of  system  performance. 

6.  Repeat  step  five  until  the  measure  of  system  performance 
stabilizes . " 


11 


Thus  the  distinguishing  feature  of  the  Monte  Carlo  method  is  the  re¬ 
petitive  execution  of  an  established  experiment  or  simulation  involving 
randomness . 

While  electronic  digital  computers  themselves  are  not  necessary  for 
the  execution  of  simulations,  they  do  offer  tremendous  speed  and  con¬ 
sistency  of  conditions  for  such  models.  Thus  the  computer  is  ideally 
suited  to  perform  the  large  number  of  repetitions  required  by  Monte 
Carlo  but  the  matter  of  randomness  presents  a  problem.  For  the  Monte 
Carlo  technique  described  above  the  necessity  of  a  truly  random  number 
is  essential.  However  Kovach  |jtef.  ej  notes  that: 

"Strictly  speaking,  the  random  number  exists  only  as  the 
result  of  a  random  process." 

While  computers,  to  include  the  TI  59,  do  possess  the  capacity  to  con¬ 
tinuously  generate  random  numbers  as  they  are  needed,  these  numbers  are 
subject  to  the  limitations  of  the  computer  and  are  not  truly  the  result 
of  a  random  process  and  hence  are  often  described  as  pseudo-random. 


a.  TI  59  RANDOM  NUMBER  GENERATOR 

R.F.  Barton  |~Ref.  describes  simulation  as  follows: 

"Simulation  is  simply  the  dynamic  execution  or  manipu¬ 
lation  of  a  model  of  an  object  system  for  some  purpose . 

Simulation  is  a  case -by-case  method  for  studying 
object  systems.  Each  case  might  be  either  a  single  trial 
or  an  entire  run.  In  either  view,  outputs  may  differ 
trial  to  trial  and  run  to  run." 

The  object  system  is  that  system  under  study  in  the  simulation. 

The  TI  59  in  its  capacity  as  an  electronic  computer  provides  the 
user  with  the  means  of  developing  and  executing  stochastic  and  nonstochastic 
simulations . 
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Barton  continues  j^ef.  ^ 


"  A  stochastic  simulation  is  one  in  which  differing  outputs 
trial  to  trial  can  be  obtained  without  changing  the  inputs 
(ignoring  random  numbers  as  inputs).  Specifically,  this 
means  that  identical  parameters,  starting  conditions,  and  in¬ 
put  time  path  values  produce  varying  outputs  trial  to  trial 
and  run  to  run. 


A  nonstochastic  simulation  is  one  in  which  the  inputs  or 
the  model  must  be  changed  to  obtain  changed  outputs.  This 
means  that  identical  model  operations,  parameters,  starting 
conditions,  and  input  time  path  values  will  produce  identical 
outputs  run  to  run.  " 

There  are  inputs  common  to  both  of  these  simulation  types.  However, 
as  alluded  to  above,  there  are  also  special  inputs  that  are  needed  to 
represent  the  chance  processes  or  stochastic  events  found  only  within 
a  stochastic  simulation.  These  special  inputs  are  random  numbers. 

The  characteristic  of  random  numbers  that  makes  them  different  from 
all  other  numbers  is  the  fact  that  the  knowledge  of  any  future  random 
number  cannot  be  enhanced  by  the  knowledge  of  any  past,  present,  or 
other  future  random  number. 


The  TX  59  with  its  master  library  module  solid  state  software  pro¬ 


gram  ML-15 ,  a  random  number  generator,  can  generate  sequences  of  uniformly 


or  normally  distributed  random  numbers  independent  of  a  simulation  pro¬ 
gram  or  within  such  a  program. 

Kovach  states  [Ref.  €~1  : 

L  _1 


" (Random)  numbers  generated  by  the  computer  are  sometimes- 
called  pseudo-random  because  they  are  subject  to  the 
limitations  of  the  computer.  Xn  a  list  of  truly  random 
numbers,  for  example,  one  would  expect  to  find  numbers 
containing  more  digits  than  can  be  obtained  in  a  computer." 


Random  numbers  produced  by  the  TI  59  ML-15  program  are  generated  by 


a  mathematical  formula.  Given  an  initial  seed  number  by  the  user,  this 
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program  will  always  produce  the  same  list  of  pseudo-random  numbers. 

Thus  if  repeatedly  initialized  with  the  same  seed  number  the  forth¬ 
coming  random  numbers  would  be  known  and  randomness  would  not  exist. 

That  is,  every  future  random  number  could  be  predicted.  Hence,  the 
randomness  of  the  numbers  produced  by  the  ML-15  program  are  as  dependent 
upon  the  user  as  the  mathematical  formula  of  the  program  itself.  It  is 
therefore  incumbent  upon  the  routine  user  of  the  ML-15  program  to  vary 
the  seed  number  used  within  denoted  limits  to  insure  genuine  pseudo¬ 
random  numbers. 

The  TI  59  ML-15  random  number  generator  program  is  listed  in 
Table  2-1.  User  instructions  for  the  ML-15  program  |Ref.  llj  are 
listed  in  Table  2-2.  Data  register  contents  are  listed  in  Table  2-3. 

While  the  program  does  provide  the  option  of  generating  uniform 
numbers  for  ranges  other  than  0-1  and  also  provides  statistical  data 
for  the  random  numbers  generated,  only  the  generation  of  the  uniform, 
range  0-1,  and  the  normal  random  numbers  will  be  discussed  further. 

1.  Uniform  Random  Number  Generator 

A  flowchart  of  the  uniform  random  number  generator,  range  0-1, 
is  displayed  in  figure  2-1.  Program  steps  000  through  054  contain  the 
following  mathematical  formula,  called  the  linear  congruential 
method  JlRef.  llj  ,  for  the  generation  of  these  numbers.  (Throughout 
this  thesis  an  asterisk  is  used  to  indicate  multiplication.) 

((24298  *  SEED  +  99991)  -u_  199017  STO  07) 

A  work  value  is  the  result  of  the  above  operation.  This  result  remains 
in  the  display  register.  The  value  199017  is  stored  in  data  register 
07,  an  ML-15  work  register. 
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TI  59  Uniform  Random  Number  Generator 


FIGURE  2-1 
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Calculation  continues. 


(INV  I NT  *  RCL  07)  STO  09 

The  integer  portion  of  the  number  resulting  from  the  previous  operation 
is  discarded,  then  the  remaining  fractional  portion  is  multiplied  by 
199017  which  was  stored  in  data  register  07.  This  product  then  is 
stored  in  data  register  09  and  becomes  the  seed  for  the  next  random 
number  calculation.  Calculation  continues. 

((RCL  09  -j-  RCL  07)  *  5  INV  2nd  log) 

Now  the  new  seed  is  divided  by  the  number  199017  which  was  stored  in 
data  register  07  during  the  first  operation.  This  quotient  is  then 
multiplied  by  the  common  antilogarithm  of  5  to  complete  the  step. 
Calculation  continues. 

(INT  -r-  5  INV  log) 

The  fractional  portion  of  the  previous  numerical  operation  is  discarded 

and  the  result  divided  by  the  common  antilogarithm  of  5.  This  quotient 

is  then  displayed  as  the  uniform  random  number,  range  0-1. 

2 .  Normal  Random  Number  Generator 

A  flowchart  of  the  normal  random  number  generator  is  displayed 

in  figure  2-2.  Program  steps  069  through  135  contain  the  following 

mathematical  formula,  called  the  direct  method  £ef.  l]  for  the 

generation  of  these  numbers. 

(seed  STO  09) 

(mean  STO  10) 

(standard  deviation  STO  11) 

Following  program  initialization,  three  data  values  are  entered  and 
stored  in  data  registers  09,  10,  11.  The  seed  number,  the  desired 
normal  distribution  mean  (mu)  and  the  desired  normal  standard  deviation 
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(sigma)  are  stored  in  these  registers  as  listed  respectively.  The  seed 
value  is  limited  as  noted  in  the  TI  user  instructions,  table  2-2. 

(SBR  DMS  STO  08  SBR  DMS) 

Initially  the  program  calls  the  previously  discussed  uniform  random 
number  generator  to  produce  such  a  number  and  then  stores  it  in  data 
register  08.  The  uniform  generator  is  called  again  to  produce  another 
uniform  random  number  which  remains  in  the  display  register  for  manipu¬ 
lation  and  is  denoted  RN.  Calculation  continues. 

(CRN  *  2  *  TT  )  COS)  *  (  (RCL  08  LnX)  *  (-2) 

The  uniform  random  number  in  the  display  is  next  multiplied  by  two  pi. 

The  cosine  of  this  product  is  then  calculated.  The  resulting  value  is 
multiplied  by  the  product  of  the  natural  logarithm  of  the  first  uniform 
random  number  (generated  early  and  stored  in  data  register  09)  multiplied 
by  negative  two.  This  product  remains  in  the  display  register  for 
manipulation.  Calculation  continues. 

("V - *  RCL  11)  +  RCL  10 

The  squareroot  of  the  previous  operation  end  value  is  multiplied  by  the 
desired  standard  deviation.  Finally,  this  product  is  added  to  the 
desired  mean,  resulting  in  the  generated  normal  random  number.  Seed 
manipulation  for  the  generation  of  successive  normal  random  numbers  is 
completed  during  the  SBR  D.MS  portion  of  the  normal  generation  program. 

3.  SBR  2nd  D.MS 

The  TI  59  ML  15  program  does  compile  statistical  data  to  allow 
computation  of  the  mean  and  standard  deviation  of  the  pseudo-random 
numbers  generated  when  using  the  normal  distribution  routine  and  the 
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uniform  distribution  routine  over  ranges  other  than  0-1.  However,  these 
data  are  not  compiled  when  SBR  2nd  D.MS  sequence  is  executed  to  pro¬ 
duce  uniformly  distributed  numbers  over  the  0-1  range.  Data  registers 
one  through  eleven  are  used  by  the  ML-15  program  to  compile  and  compute 
these  statistics.  Hence,  if  this  program  is  called  to  produce  normal 
random  numbers  within  a  larger  program,  such  as  a  simulation,  the  use 
of  these  eleven  data  registers  must  be  forgone.  Yet,  if  the  ML-15  pro¬ 
gram  is  called  only  to  produce  uniform  random  numbers  over  the  range 
0-1,  only  data  registers  seven  and  nine  are  used  by  ML-15,  freeing 
nine  registers  for  other  use .  This  aspect  of  the  TI  59  ML-15  program 
must  be  carefully  considered  when  utilizing  it  is  as  a  subroutine  within 
another  program. 

C .  COMBAT  MODELS 

Today  there  are  considered  to  be  three  types  of  combat  models  in 
use;  war  games,  pure  simulations  and  analytical  models.  War  games  are 
models  and  games  2]  ■ 

"...  in  which  individuals  simulating  decision  makers 
in  real  life  use  their  judgement  to  perform  the 
decision  functions  in  the  model." 

A  war  game  may  include  automation  to  assist  in  the  processing  of  data 
and  the  generation  of  random  numbers  to  determine  the  outcome  of  certain 
chance  events.  A  war  game  may  also  be  a  player-assisted  simulation 
where  players  provide  input  to  a  computer  model  based  on  output  (read¬ 
outs)  during  a  simulated  battle.  In  comparison  with  the  other  models, 
war  games  appear  to  be  more  realistic,  involve  greater  player  interac¬ 
tion,  are  less  automated,  require  much  more  time  to  run,  more  resources 
and  involve  a  smaller  degree  of  abstraction. 
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Simulation  combat  models  are  models 


[*ef.  2J  : 


"...  which  run  completely  without  human  intervention. 

In  this  type  of  model  events  in  the  different  combat 
processes  are  based  on  predetermined  rules  which  are 
programmed  into  the  automated  evaluation  procedure." 

Combat  models  of  this  type  generally  contain  a  significant  number  of  the 

important  stochastic  elements  of  combat  in  an  attempt  to  simulate  real 

battle.  These  models  use  probability  distributions  for  the  many  chance 

input  variables  and  produce  probability  distributions  as  results. 

They  utilize  the  Monte  Carlo  technique ,  repeatedly  sampling  all  input 

distributions  in  the  programmed  sequence  to  produce  a  distribution  of 

probable  battle  results  for  each  set  of  input  data. 


Analytical  models  are  models  |~Ref.  2J  : 


"...  comprised  of  sets  of  mathematical  equations  as 
models  of  all  the  basic  events  and  activities  in  the 
process  being  described  and  an  overall  assumed 
mathematical  structure  of  the  process  into  which,  the 
event  or  activity  descriptions  are  integrated.  " 

While  analytical  models  are  the  most  time  efficient  they  are  also  the 
most  abstract  and  difficult  to  understand.  As  with  the  pure  simulation, 
there  is  no  human  intervention  when  an  analytical  model  is  used. 

All  three  models  represent  abstractions  of  the  real  world.  The 
models  themselves  can  be  observed  more  conveniently  than  the  real  world 
and  theories  about  the  real  world  can  be  developed  by  studying  the  re¬ 
sults  of  these  models.  Subsequently,  these  theories  can  be  used  to 
make  predictions  about  real  world  events. 

Each  model  type  has  strengths  and  limitations ,  some  noted  above  and 


others  listed  in  table  2-4 


Ref.  2]  . 
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1.  Pure  Combat  Simulations 


Pure  combat  simulations  are  normally  viewed  as  production  tools , 
using  Monte  Carlo  techniques  to  obtain  results  enabling  the  prediction 
of  future  system  performance.  But  because  the  real  world  is  so  complex 
and  interactive,  attempts  to  model  every  detail  of  a  large  system  in  a 
pure  simulation  and  to  include  every  element  that  may  influence  the 
system  can  result  in  simulations  so  large  and  so  complex  that  they 
are  understood  only  by  their  developers  and  not  by  other  users  or 
decision  makers.  To  avoid  this  complex  dilemma,  analytical  models 
can  be  used  to  represent  elements  of  the  system  being  modeled  instead 
of  simulating  the  element  itself  and  its  inherent  stochastic  processes 
with  every  trial.  This  technique  has  been  followed  to  some  degree  in 
the  pure  simulation  model  presented  in  Chapter  III  where  the  calculation 
of  detection  probability  is  an  analytical  model  with  results  based 
on  target  range. 

It  should  be  noted  that  few,  if  any,  simulation  models  ever  com¬ 
pletely  include  all  those  elements  and  events  that  affect  the  system (s) 
under  study.  Reference  2  points  out  that: 

".  .  .a  model  is  always  incomplete,  with  only  those 
aspects  represented  that  we  believe  we  know  well 
enough  to  model  and  that  we  consider  important  in 
the  issues  to  be  examined  with  the  model.  Obviously, 
models  tend  to  be  as  simple  and  concise  as  our 
knowledge  of  the  activity  warrants . " 

This  is  reflected  in  the  model  presented  in  Chapter  III.  While  all  the 

factors  affecting  system  performance  have  not  been  directly  simulated 

they  are  included  either  as  analytical  models  or  as  given  in  the 

scenario. 
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Finally,  a  point  to  be  stressed  is  that  simulations  need  not  be 
large  to  be  useful,  nor  require  the  use  of  a  large  electronic  digital 
computer  to  be  credible .  Using  large  computers  for  large  problems  and 
small  computers  for  small  problems  is  a  rule  of  thumb  that  may  over¬ 
state  the  case  but  certainly  does  not  exaggerate  it.  Use  of  the  TI  59 
as  a  computer  to  tackle  the  problem  set  forth  in  Chapter  III  is  an 
example  of  matching  the  problem  to  the  computing  power  required.  It 
is  also  an  excellent  example  of  the  computing  power  of  the  TI  59. 
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TI  59  PROGRAM  ML-15 


USER  INSTRUCTIONS 
(MASTER  LIBRARY  MODULE) 


STEP 

PROCEDURE 

ENTER 

PRESS 

DISPLAY 

1 

SELECT  PROGRAM 

2d 

PGM  15 

T 

INITIALIZE 

2d 

E' 

0 

3 

ENTER  RANDOM  NUMBER  SEED 

SEED 

E 

SEED 

(0  S  SEED  £  199017) 

FOR  NORMAL  DISTRIBUTION 

4 

ENTER  DESIRED  MEAN 

s 

A 

5 

ENTER  DESIRED  STANDARD 

DEVIATION 

a- 

B 

6 

GENERATE  RANDOM  NUMBER 

(REPEAT  AS  NEEDED) 

2d 

c1 

RANDOM 

NUMBER 

FOR  UNIFORM  DISTRIBUTION 

RANGE  (0,  1) 

7 

GENERATE  RANDOM  NUMBER 

(REPEAT  AS  NEEDED) 

SBR 

2d  D .MS 

RANDOM 

NUMBER 

TABLE  2-2 
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TI  59  PROGRAM  ML-15 


DATA  REGISTER  CONTENTS 


REGISTERS: 

R 

00 

R0l 

ly 

R._ 

02 

zy2 

R03 

N 

R04 

lx2 

R05 

lx 

R06 

Z  xy 

R07 

WORK  REGISTER 

R08 

WORK  REGISTER 

**09 

SEED 

R10 

A,  % 

*11 

B,«r 

TABLE  2-3 
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A 


War  Games 


Simulations 


Intelligent  play  of  decision  maker. 

Intelligent  use  of  intelligence. 

Can  plan  engagement  and  moves  in 
advance . 

Adapts  maneuver  to  situation. 


Can  play  many  tactical  situations 
(employment,  penetration,  etc.). 

Insight  is  gained  by  understanding 
the  particular  rationale  used  in 
the  decision  process  in  single 
situations  analyzed. 

Controller  determines  existence 
of  engagement  and  pace  of  play. 

Very  slow  and  costly  in  resources. 

Very  few  situations  can  be 
examined. 


Stylized  decision  routines, 
usually  fixed  throughout  game . 

Very  limited  use  of  intelligence. 

Very  limited  planning  horizon. 

Very  limited  adapted  maneuver 
routines . 

Usually  stylized  maneuver,  limited 
change  in  formations. 

Insight  is  gained  by  repeating 
the  analysis  in  many  situations 
using  different  values  for  key 
parameters . 

Predetermined  scenarios  and 
engagement  rules  —  combat  very 
intense . 

Faster  to  run  after  completely 
developed. 

Many  situations  can  be  played 
and  the  sensitivity  of  key 
variables  can  be  tested. 


Greatest  visibility  for  the  user. 

Can  include  direct  involvement  by 
user . 


^Reference 


1,  FM  44-90,  Headquarters  Department  of  the  Army,  1977 
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III.  combat  simulation  USING  TI  59 


This  is  a  probabilistic  duel  simulation  model,  a  pure  simulation  of 
a  combat  air  battle,  designed  to  reflect  the  characteristics  of  the 
Improved  HAWK  air  defense  artillery  system  (battery)  in  the  manual 
mode  under  attack  by  a  single  cruise  missile  of  sustained  altitude, 
speed  and  direction. 

The  program  scenario  and  engagement  rules  are  predetermined  with 
no  user  input  once  the  simulation  run  has  begun.  Insight  may  be  gained 
and  the  sensitivity  of  key  variables  tested  dependent  upon  the  use  of 
different  values  entered  by  the  user  for  these  variables  during  program 
initialization.  (Variables  listed  under  E  below.) 

The  model  provides  IHAWK  system  status ,  target  engagement  events  and 
battle  results  as  they  are  determined/ occur .  Only  two  battle  results 
are  possible:  a  "KILL"  of  the  cruise  missile  or  a  unit  "PENETRATION" 
by  the  cruise  missile. 

Given  the  operational  ready  rates  of  the  major  subsystems  of  the 
IHAWK  system,  the  Pssk  (probability  of  single  shot  kill)  and  the  mean 
and  standard  deviation  of  lock-on-to-target  times,  the  model  samples 
from  the  uniform  and  normal  distribution  to  determine  system  status, 
IHAWK  missile  kill  or  no  kill  and  lock-on  times.  Target  detection  is 
modeled  as  a  function  of  target  range  and  is  represented  as  a  linear 
relationship  in  the  simulation. 

This  simulation  was  developed  to  exhibit  the  computing  power  of  the 
TI  59  and  to  determine  if  one  parameter  under  the  control  of  the  IHAWK 
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battery  commander  could  significantly  affect  air  battle  results.  This 
parameter  was  the  tracking  radar  " lock-on- to-target "  time  which  is  a 
function  of  operator  training  given  (1)  a  manual  mode  operation,  and 

(2)  perfect  equipment.  The  sensitivity  of  battle  results  to  varying 
lock-on  times  is  listed  under  F  below. 

A.  I HAWK  SYSTEM 

To  be  effective,  an  air  defense  system  must  be  able  to  detect, 
identify,  engage  and  destroy  an  airborne  target.  The  I HAWK  system  can 
engage  and  destroy  a  full  spectrum  of  threat  aircraft  and  missiles  operat¬ 
ing  throughout  a  wide  range  of  tactical  speeds  and  altitudes.  It  can 
engage  a  multiple  target  threat  as  well  as  single  targets.  The  system 
is  effective  from  ground  level  up  to  altitudes  of  about  48000  feet  and 
out  to  ranges  of  about  40  kilometers.  The  system  can  operate  at  night, 
under  all  conditions  of  weather  and  reduced  visibility.  It  can  function 
effectively  in  an  ECM  (electronic  countermeasures )  environment  and  is 
mobile  using  organic  unit  vehicles  or  helicopters  [■Ref.  lj  .  However, 
some  adverse  weather  and  heavy  ECM  may  diminish  some  system  capabilities. 

This  simulation  does  not  model:  (1)  weather 

(2)  detection  ECM 

(3)  visibility 

(4)  system  mobility 

1.  Detection 

Target  detection  is  accomplished  by  either  the  improved  pulse 
acquisition  radar  (IPAR)  or  the  improved  continuous  wave  acquisition 
radar  (ICWAR) ,  or  both.  The  IPAR  can  detect  low  to  medium  altitude 
targets  out  to  ranges  in  excess  of  100  kilometers  while  the  ICWAR  can 
detect  targets  at  very  low  altitudes  with  ranges  in  excess  of  60 


31 


1C  WAR 


kilometers.  Operating  on  the  doppler  principle,  the  ICWAR  sees  only 
very  low  moving  objects  and  not  stationary  objects  on  the  ground. 
Detection  may  be  reduced  if  terrain  features  such  as  hills,  trees  and 
buildings  mask  the  radar's  view  of  the  approaching  target.  Even  with 
completely  level  terrain,  the  earth's  curvature  causes  a  reduction  in 
the  detection  and  tracking  ranges  of  the  system  against  very  low-altitude 
targets.  Also,  evasive  maneuvers  by  threat  targets  may  reduce  the  de¬ 
tection  and  tracking  ranges  and  increase  system  response  time,  thereby 
reducing  the  effective  or  intercept  range.  This  simulation  does  not 

model:  (1)  terrain  features 

(2)  curvature  of  the  earth 

(3)  evasive  action  by  targets 

(4)  pulse  detection  or  continuous  wave  detection  per  se 
This  model  assumes: 

(1)  clear  weather 

(2)  no  detection  electronic  countermeasures 

(3)  line  of  sight  (LOS)  exists  between  radar  and  target 

(4)  flat  desert  terrain 

( 5 )  nonmaneuvering  target 

(6)  only  one  attacking  target  exists 

(7)  target  is  a  cruise  missile  of  constant  speed 
and  constant  altitude 

(8)  detection  is  a  function  of  target  range 

The  probability  of  detection  is  modeled  during  each  sweep  of  the  radar 
as  a  linear  function  of  target  range  from  the  battery  as  follows: 

For  the  IPAR  the  P^et  =  ^-•25  4*  65)  *  Target  RangeJ  +  1.0 

For  the  ICWAR  the  PDET  =  £(-.5  -i—  65)  *  Target  RangeJ  +  1.0 

The  probability  of  detection  is  calculated  every  three  seconds  of 
simulated  time.  This  is  based  on  the  radar  rotation  rate  of  20 
revolutions  per  minute.  That  is,  every  three  seconds  each  radar 
takes  a  360  degree  glimpse  of  the  horizon.  The  radars  are  slaved  to 
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each  other  and  rotate  in  synchronization.  Additionally,  the  IPAR  is 
modeled  to  detect  only  targets  from  5000  to  40000  feet  in  altitude  while 
the  ICWAR  detects  targets  from  1  to  8000  feet  in  altitude.  Thus  the 
battery's  very  low  and  low  to  medium  detection  capability  is  dependent 
upon  the  operational  status  of  these  radars  as  noted.  Targets  above 
40,000  feet  cannot  be  detected  in  this  simulation.  The  operational 
ready  rates  of  these  two  radars  has  been  arbitrarily  set  at  .65  (ICWAR) 
and  .95  (IPAR). 

2 .  Identification 

Identification  of  any  potential  target  is  accomplished  by  means 
of  the  identification,  friend  or  foe  (IFF),  equipment  of  the  I HAWK 
system  and/or  other  established  hostile  criteria.  If  the  target  can¬ 
not  be  positively  identified  in  this  simulation  because  of  a  non- 
operational  IFF,  the  target  speed  and  altitude  is  checked  to  determine 
target  status  (foe  or  not  foe) .  That  is,  if  the  IFF  is  nonope rational 
and  if  the  target  is  below  5000  feet  altitude  and  greater  than  550  KMPH 

in  speed,  it  is  identified  as  a  foe;  otherwise,  it  is  not  a  foe. 

This  model  assumes;  (1)  once  identified  as  a  friend,  always  a  friend, 

(2)  once  identified  as  a  foe,  always  a  foe. 

The  operational  ready  rate  of  the  IFF  has  been  arbitrarily  set  at 

95  percent.  This  model  does  play  IFF  accuracy  to  the  degree  that  an 

operational  IFF  will  be  in  error  two  percent  of  the  time.  That  is,  a 

foe  will  be  shown  to  be  a  friend  two  percent  of  the  time.  This  model 

assumes  operational  IFF  accuracy  to  be  98  percent. 
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o  Identify 


Figure  3-2 


-  -  IFF  CHALLENGE  — - - - 

_  CORRECT  RESPONSE  IF"'*"0* 

NO  RESPONSE  OR  .NCORRECT  RESPONSE! 
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3 .  Engagement/Destruction 

The  IHAWK  battery  has  two  firing  sections,  each  of  which 
contains  a  target  tracking  radar  which  is  called  an  improved  high  powered 
illuminating  radar  (IHIPIR) ,  and  three  launchers  with  three  missiles 
each.  If  a  potential  target  is  determined  to  be  a  foe,  it  is  assigned 
to  one  or  both  of  these  firing  sections.  The  tracking  radar  of  these 
sections,  under  the  control  of  a  fire  control  operator  (enlisted 
personnel) ,  attempts  to  lock-on  to  the  approaching  hostile  target.  The 
operator  directs  the  automatic  box  search  of  the  radar  in  the  azimuth  and 
expected  elevation  of  the  target  in  attempting  this  target  lock.  The 
operator's  ability  to  achieve  a  target  lock  is  a  function  of  his  train¬ 
ing  given  the  condition  mentioned  above.  The  time  that  elapses  in 
attempting  this  target  lock  is  extremely  important.  Target  engagement 
cannot  continue  without  target  lock  and  the  longer  it  takes  to  achieve 
target  lock  the  closer  the  target  moves  toward  the  battery,  reducing  the 
intercept  range.  Tracking  radar  lock-on  times  are  assumed  to  be  normally 
distributed.  After  target  lock  has  been  achieved  and  the  target  is  in 
range,  one  or  two  missiles  are  fired  on  order  from  the  unit  tactical 
control  officer  in  a  battery  control  van.  Engagement  is  continued 
until  the  target  is  destroyed  or  until  engagement  is  no  longer  possible. 
This  simulation  models  each  firing  section  as  an  entity.  After  a 
target  has  been  declared  a  foe  and  assigned  to  one  or  both  firing 
sections,  this  model  simulates  the  target  lock-on  time  by  utilizing  the 
TI  59  normal  random  number  generator  based  on  a  normal  mu  and  sigma 
input  by  the  user  during  initialization. 
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TIME 


HAWK 

SITE 


RANGE 


FLIGHT  TIME 


INTERCEPT 


FIRE 


“  TIME  I 

TRACK  DETECTION 


The  range  at  intercept  is  determined  by  the  range  at  which  detection  and  tracking 
(lock-on  by  the  HIPIR)  occur  and  on  system  response  time 


In  this  model,  targets  are  engaged  that  are: 

(1)  declared  to  be  foe 

(2)  less  than  40  KM  from  the  battery 

(3)  greater  than  8  KM  from  the  battery 

This  model  assumes: 

(1)  two  independent  firing  sections 

(2)  salvo  fire  occurs  if  both  sections  are  operational  and 
shoot-look-shoot  if  only  one  section  is  operational 

(3)  firing  continues  until  kill  or  penetration 

(4)  penetration  means  that  the  target  is  8  KM  or  less  from 
the  battery 

(5)  lock-on-to-target  time  is  a  function  of  operator 
training  and  is  normally  distributed 

(6)  Pgs!c  =  -75  (arbitrarily  set) 

The  operational  ready  rate  of  each  firing  section  has  been  arbitrarily 
set  at  75  percent. 

4.  Target 

The  target  for  this  model  is  assumed  to  be  a  hostile  cruise 
missile  that  flies  straight  in  toward  the  battery  at  a  constant  speed 
and  altitude  as  established  by  the  user  during  initialization.  The 
initial  range  of  the  target  is  also  a  user  input.  The  lethality  of 
the  missile  warhead  is  assumed  to  be  such  that  any  successful  penetra¬ 
tion  by  the  missile  to  within  8  KM  or  less  of  the  battery  before 
destruction  is  considered  a  total  penetration  of  the  battery  defended 
area.  Therefore,  the  target  must  be  destroyed  before  8  KM  to  score  a 
kill.  Additionally,  a  target  will  not  be  engaged  after  detection  until 
it  is  less  than  40  KM  from  the  battery  and  no  further  missiles  will  be 
fired  at  the  target  once  it  is  within  8  KM  of  the  battery.  The  target 
speed  has  a  lower  bound  of  100  KMPH  but  no  upper  bound.  Only  targets 
between  1  and  40000  feet  in  altitude  can  be  detected  and  are  thus  the 


altitude  bounds.  Finally,  all  targets  are  hostile  and  will  be  engaged 


unless  erroneously  identified  as  friendly  or  not  foe,  resulting  in  a 
free  penetration. 

5 .  Time 

This  simulation  is  a  time  step  model ,  updating  all  battery 
events  and  functions  every  three  seconds  of  simulated  time.  This  three- 
second  interval  stems  from  the  rotation  rate  of  the  detection  radars, 

20  revolutions  per  minute  or  one  complete  rotation  (scan  of  the  horizon) 
every  three  seconds. 

B.  MACRO  FLOW  CHART 

The  enclosed  macro  flow  chart,  figure  3-4,  depicts  the  general  flow 
of  the  simulation  logic  from  start  to  either  penetration  or  kill. 

First  the  model  determines  if  a  detection  capability  exists.  This 
could  be  one  or  both  of  the  detection  radars.  Using  the  internal  random 
number  generator  of  the  TI  59  for  a  0-1  uniform  distribution,  two 
random  numbers  are  drawn  and  compared  with  the  detection  radar  opera¬ 
tional  ready  rates.  If  the  random  number  is  less  than  the  rate,  the 
radar  is  operational;  otherwise,  it  is  nonoperational.  If  no  detection 
capability  exists  the  simulation  is  terminated  by  a  penetration  of  the 
defended  area  by  the  target. 

Detection  of  the  target  is  based  partly  on  the  formulas  set  forth 
above  and  results  of  the  0-1  range  uniform  random  number  generator. 

The  probability  of  detection  is  based  on  the  range  of  the  approaching 
target  and  is  recalculated  every  three  seconds.  The  probability  of 
detection  for  each  radar  is  compared  with  a  generated  random  number 
between  0  and  1.  If  the  random  number  is  less  than  the  probability 
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the  target  is  detected;  otherwise,  it  is  not.  As  the  target  moves 
closer  to  the  unit,  the  probability  of  detection  increases  linearly. 
However,  detection  is  not  based  on  range  alone.  The  target  must  be 
within  the  detection  altitude  capability  of  the  operational  system. 

For  instance,  if  the  battery's  detection  capability  rests  solely  on  the 
IPAR  because  of  a  nonoperational  ICWAR  and  if  the  approaching  target  is 
at  an  altitude  of  less  than  5000  feet,  detection  will  never  occur  and 
a  penetration  will  result. 

Identification  of  the  target  is  determined  by  either  the  IFF  or 
a  combination  speed  and  altitude  envelope  if  the  IFF  is  nonoperational. 
Operational  status  of  the  IFF  is  determined  by  comparing,  again,  a 
generated  random  number  from  a  0-1  range  uniform  distribution  against 
the  operational  ready  rate  of  the  IFF.  If  the  random  number  is  less 
than  the  ready  rate,  the  IFF  is  operational;  otherwise,  it  is  non¬ 
operational.  If  the  IFF  is  nonoperational  and  if  the  unit  has  a 
tracking  capability,  the  target  speed  and  altitude  can  be  checked 
against  an  established  hostile  criteria.  If  the  target  is  found  to  be 
below  5000  feet  in  elevation  and  above  550  KMPH  in  speed,  it  is 
designated  a  foe;  otherwise,  it  is  considered  to  not  be  a  foe.  Only 
targets  positively  identified  as  foes  are  engaged. 

IFF  positive  identification  of  the  hostile  cruise  missile  as  a  foe 
is  determined  by  comparing  a  random  number  generated  from  a  0-1  range 
uniform  distribution  against  the  IFF  accuracy  rate  of  98  percent.  If 
the  random  number  is  less  than  the  accuracy  rate,  the  cruise  missile 
is  correctly  identified  as  a  foe;  otherwise,  it  is  erroneously  classi¬ 
fied  as  a  friend  resulting  in  a  penetration  of  the  defended  area. 
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The  battery's  engagement  capability  lies  with  its  two  firing 
sections,  each  tracking  and  firing  on  approaching  hostile  targets. 

Again,  two  random  numbers  drawn  from  the  0-1  uniform  distribution  are 
compared  to  the  operational  ready  rates  of  the  firing  sections  to  de¬ 
termine  the  status  as  in  the  previous  subsystems  examples.  If  an 
engagement  capability  is  determined  to  exist,  the  process  continues  on 
to  direct  missile  firings  at  the  target  following  target  lock  by  the 
tracking  radar (s) .  As  mentioned  earlier,  this  lock-on  time  is  a 
function  of  operator  training  in  the  manual  mode  and  is  normally 
distributed.  For  highly  trained  operators  the  mean  is  assumed  to  be 
ten  seconds  with  a  standard  deviation  of  five  seconds.  The  lock-on 
time  for  each  section  is  determined  by  two  random  numbers  generated  by 
the  TI  59  random  number  generator  from  a  normal  distribution  with  mean 
and  standard  deviation  determined  by  the  user  during  program  initializa¬ 
tion.  (Any  random  number  less  than  zero  is  discarded  and  another 
generated  to  avoid  negative  times.)  The  length  of  the  lock-on  period 
directly  affects  the  resultant  target  intercept  range.  If  the  lock-on 
time  is  sufficiently  long  the  missile  firing  is  delayed  and  the  proba¬ 
bility  of  a  penetration  is  likely.  (The  determination  of  air  battle 
results  to  varying  lock-on  times  provided  the  basis  for  the  development 
of  this  simulation,  though  other  variables  of  the  model  can  easily  be 
tested  for  outcome  sensitivity.)  Again,  targets  are  not  engaged  until 
less  than  40  KM  from  the  battery  and  no  missiles  are  fired  after  the 
target  is  3  KM  or  less  from  the  battery. 

Missile  effectiveness  after  firing  is  determined  when  the  IHAVJK 
missile  range  equals  or  exceeds  the  cruise  missile  range  from  the 
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battery.  A  random  number  from  a  0-1  range  uniform  distribution  is 
generated  for  each  missile  fired  and  checked  against  the  Pg^. 

Random  numbers  less  than  the  Pgg^  result  in  kills  while  all  others 
result  in  no  kills.  Engagement  of  the  target  continues  until  a  kill 
or  penetration  is  registered. 

C.  MICRO  FLOWCHART 

The  micro  flowchart  in  appendix  A  depicts  the  detailed  flow  of 
processing  throughout  the  simulation  from  start  to  kill  or  penetration. 

D.  SUBROUTINES,  LABELS,  FLAGS,  DATA  REGISTERS  AND  PROGRAM  MEMORY  STEPS 

1.  Subroutines  and  Labels 

This  simulation  uses  49  of  the  72  labels  available  for  pro¬ 
gramming  on  the  TI  59.  Of  these  49,  14  are  subroutines.  The  remaining 
labels  are  used  to  identify  sections  of  the  program  and  to  direct 
action  to  these  sections  during  simulation  runs.  A  complete  listing 
of  all  labels  is  displayed  in  table  3-1  with  the  subroutines  marked 
by  an  asterisk.  Comments  on  selected  labels  and  a  printout  of  the 
entire  program  is  enclosed  in  appendix  A. 

2 .  Flags 

Nine  of  ten  available  flags  are  used.  As  the  I HAWK  equipment 
and  target  friend/foe  status  is  determined,  this  model  uses  TI  59  flags 
to  maintain  a  record  of  the  system  and  target  status.  These  flags 
subsequently  direct  the  flow  of  processing  and  determine  actions  to 
be  taken  within  the  simulation  engagement. 

Flag  1  set  means  the  target  has  been  detected.  Flag  2  set 
means  the  target  is  a  friend.  Flag  3  set  means  the  ICWAR  and  the  IPAR 
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are  operational.  Flag  4  set  means  only  the  ICWAR  is  operational. 

Flag  5  set  means  only  the  IPAR  is  operational.  Flag  6  set  means  the 
IFF  is  operational.  Flag  7  set  means  that  alfa  firing  section  is 
operational  and  bravo  firing  section  is  nonoperational .  Flag  8  set 
means  that  bravo  firing  section  is  operational  and  alfa  firing  section 
is  nonoperational.  Flag  0  set  means  that  both  firing  sections  are 
operational . 

3 .  Data  Registers 

The  TI  59  memory  storage  area  is  initially  partitioned  to  provide 
60  data  storage  registers  and  480  program  storage  locations.  However, 
the  user  can  repartition  the  memory  storage  area  to  suit  his  particular 
programming  needs.  The  IHAWK  simulation  requires  exactly  800  program 
memory  locations  and  20  data  storage  registers.  Within  the  TI  59  there 
are  a  total  of  120  registers  to  be  used  for  data  storage  and  program 
locations.  While  each  register  can  store  only  one  datum  point,  each  can 
store  eight  program  instructions  or  steps.  Thus  3  *  60  =  480  program 
locations  which  are  initially  available  as  mentioned  above.  Reparti¬ 
tioning  the  core  120  registers  is  done  in  increments  of  ten.  Hence,  to 
get  the  800  program  steps  for  the  IHAWK  simulation  100  core  registers 
are  needed.  This  leaves  exactly  the  20  needed  for  data  storage. 

To  partition  the  storage  area,  the  number  of  sets  of  10  data 
registers  needed  is  entered  and  2nd  OP  17  pressed.  Thus  for  the  IHAWK 
simulation,  twenty  data  registers  are  available  after  the  initial  repar¬ 
titioning  by  pressing  2  2nd  OP  17.  The  registers  and  their  contents  are 
listed  on  the  following  page. 
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Roo 

**01 


02 

R03 

**94 

R05 

R06 

R07 

R08 

**09 

R10 


R12 

R13 

R14 

r15 

R16 


18 


'19 


4. 


Target  range.  Entered  by  user, 

Rq2_  through  RQg  are  used  by  the  TI  59  random  number  generator 
program. 


SEED  for  random  number  generator.  Entered  by  user. 

Mean  lock-on-to-target  time.  Entered  by  user. 

Standard  deviation  of  lock-on-to-target  time.  Entered  by  user. 

IHAWK  missile  range  from  battery.  Initially  zero. 

I HAWK  missile  range  from  battery.  Initially  zero. 

Probability  of  detection  work  register. 

Target  speed.  Entered  by  User. 

Target  altitude.  Entered  by  User. 

Target  range  work  register.  Not  entered  by  user. 

Simulation  trills  or  runs  to  be  completed.  Entered  by  user. 

Simulated  time  in  seconds  for  each  trial.  Initially  zero  for 
each  trial. 

Program  Memory  Steps 

There  are  800  program  steps  available.  All  800  program  memory 
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steps  are  used  in  this  program. 


USER  INSTRUCTIONS 


E,  . 

The  enclosed  user  instructions,  table  3-2,  provide  the  necessary- 
steps  to  initiate  a  sequence  of  simulation  runs.  The  enclosed  printout 
results,  table  3-3,  indicate  the  22  possible  print  statements  that 
may  occur  during  the  simulation.  A  sample  of  data  input  and  simulation 
run  results  are  displayed  in  tables  3-4  and  3-5. 

Two  steps  of  the  user  instructions  warrant  further  comment. 

1.  Step  2  Clear  Data  Registers. 

Instead  of  clearing  all  data  registers  the  user  may  wish  to 
clear  selected  registers  when  repeating  simulation  runs  as  in  the  case 
of  sensitivity  analysis  work.  In  this  instance  the  user  may  just  clear 
registers  Rg^,  R1?,  Ri3'  Ri7'  R18'  Rig  and  enter  the  desired  values. 

Rgg  need  not  be  reentered  as  the  program  automatically  changes  the 
seed  after  each  random  number  is  generated.  If  the  user  does  clear 
all  data  registers  with  2d  CMs  the  user  must  then  enter  an  entirely 
different  seed  in  Rgg  within  the  bounds  noted. 

2.  Step  9  Check  Data  Register  Content  \ 

This  step  is  a  quick  safeguard  for  the  user  to  ensure  that  the 
simulation  run  is  based  on  the  correct  parameter  values.  This  step 
provides  a  complete  listing  of  the  20  data  registers  with  contents  for 
review  prior  to  the  final  user  step. 

F.  SENSITIVITY  ANALYSIS  RESULTS  FOR  LOCK-ON  TIMES 

Four  hundred  simulation  runs  were  made  with  four  lock-on  mean 
values:  10,  20,  30  and  40  seconds.  In  each  case  the  lock-on  time 
standard  deviation  was  five  seconds.  For  each  simulation  run  the 


target  was  initially  set  at  a  range  of  30  :<M  in  X,.,,  target  speed  was 
1500  KMPH  and  target  altitude  was  7000  feet.  An  initial  random  number 
seed  was  entered  for  run  number  one  but  no  further  user  seeds  were 
provided,  thus  leaving  seed  manipulation  to  the  program. 

The  results  listed  below  indicate  that  air  battle  results  are  in¬ 
deed  sensitive  to  target  lock-on  times. 


LOCK-ON  TIME 

PENETRATION 

KILL 

mu  sigma 

10  5 

11% 

39% 

20  5 

14% 

36% 

30  5 

26% 

14% 

40  5 

44% 

56% 

The  results  indicate  a  significant  increase  (12  percent)  in  defended 
area  penetrations  for  a  mu  of  30  under  the  present  scenario.  This  trend 
continues  at  an  apparent  exponential  rate.  With  a  mu  of  40  seconds, 
defended  area  penetrations  increase  another  13  percent.  Eased  on  these 
results  it  appears  advisable  to  maintain  such  a  state  of  operator 
training  that  the  mean  target  lock-on  times  be  twenty  seconds  or  less 
with  as  little  deviation  among  the  operators  as  possible.  Furthermore, 
it  seems  that  for  the  extra  training  assumed  to  be  required  to  reduce 
mean  lock-on  times  from  20  to  10  seconds  there  appears  to  be  only  a 
small  marginal  reward  in  the  reduction  of  defended  area  penetrations 
(3  percent) . 

G.  RECOMMENDATIONS 

While  the  intent  throughout  the  thesis  is  to  exemplify  the  capabili¬ 
ties  of  the  TI  59  as  a  viable,  real  world  analytical  tool,  the  results 


of  the  TI  59  simulation  lend  insight  into  an  area  that  requires 


farther  investigation,  that  being  IHAWK  target  lock-on  times.  While 
only  the  lock-on  times  themselves  were  varied  for  this  simulation,  other 
important  scenario  parameters  should  be  varied  to  acquire  an  improved 
understanding  of  how  air  battle  results  can  be  affected  by  lock-on  times. 

Future  enhancements  of  these  results  would  include  a  significant 
increase  in  simulation  runs  for  a  wide  variety  of  scenario  parameter 
settings.  While  this  TI  59  model  allows  certain  parameter  variations 
during  program  initialization,  other  parameters  such  as  acquisition 
radar  altitude  detection  capabilities  can  be  varied  with  only  minor 
adjustments  to  the  program. 

Regardless  of  whether  future  simulation  studies  are  conducted  using 
this  TI  59  model  or  a  facsimile  on  another  computer,  the  results 
above  warrant  further  research  in  this  area. 
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USER  INSTRUCTIONS 


STEP 

PROCEDURE 

ENTER 

1 

REPARTITION 

2 

800  PROGRAM  MEMORY 
LOCATIONS 

20  DATA  MEMORY 

REGISTERS 

2 

Enter  magnetic  program 
cards ,  sides  1  through 

4 

3 

CLEAR  DATA  REGISTERS 

4 

RESET  ALL  FLAGS  AND 

CLEAR  ALL  SUBROUTINE 

RETURN  REGISTERS 

5 

ENTER  I HAWK  SYSTEM  DATA 

TRACKING  RADAR 
"LOCK-ON"  TIME 

V- 

6 

ENTER  CRUISE  MISSILE 

DATA 

3ANGE  (KM) 

LC<  KM  <100 

RANGE  (KM) 

10  <  KM  <  100 

SPEED 

ALTITUDE 

(FT) 

7 

ENTER  SEED  FOR  RANDOM 
NUMBER  GENERATOR 
(OSSa  199017) 

SEED 

8 

ENTER  DESIRED  NUMBER 

OF  SIMULATION  RUNS 

* 

PRESS 

DISPLAY 

2d  OP  17 

799.19 

2nd  CMs 

(NO  CHANGE) 

RST 

(NO  CHANGE) 

STO  10 

STO  11 

GTO  777  LRN 

LRN 

STO  00 

RANGE 

STO  15 

SPEED 

STO  16 

ALTITUDE 

STO  09 

SEED 

Table  3-2.1 
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STO  18 


USER  INSTRUCTIONS 

STEP 

PROCEDURE 

ENTER 

PRESS 

DISPLAY 

9 

10 

CHECK  DATA  REGISTER 

CONTENT 

START  SIMULATION  RUN 

*  PC-100  C  PRINTER  REQU] 

■ 

0 

RED 

RST 

INV  2d 

LIST 

R/S 

0 

(SEE 

RESULTS 
POSSIBLE) * 

Table  3-2.2 
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RESULT  PRINTOUTS 


NUMBER  PRINTED 
1 

3 

4 

5 

6 

7 

8 

9 

10 
11 

12 

14 

15 

17 

18 

23 

24 

25 


MEANING 


CRUISE  MISSILE  HAS  BEEN  DETECTED  AT  (RANGE  GIVEN  IN 
KMs  FROM  UNIT) 

LOW  ALTITUDE  DETECTION  RADAR  (ICWAR)  AND  MEDIUM 

ALTITUDE  DETECTION  RADAR  (IPAR)  ARE  BOTH  OPERATIONAL 

ONLY  ICWAR  IS  OPERATIONAL 

ONLY  IPAR  IS  OPERATIONAL 

IDENTIFICATION  FRIEND-OR-FOE  (IFF)  IS  OPERATIONAL 

ALFA  FIRING  SECTION  IS  OPERATIONAL,  BRAVO  FIRING 
SECTION  IS  NONOPERATIONAL 

BRAVO  FIRING  SECTION  IS  OPERATIONAL,  ALFA  FIRING 
SECTION  IS  NONOPERATIONAL 

BOTH  FIRING  SECTIONS  ARE  OPERATIONAL 

CRUISE  MISSILE  IS  IDENTIFIED  AS  A  FOE 

BRAVO  FIRING  SECTION  IS  FIRING  ONE  MISSILE  AT  A  TARGET 
<40  KM  FROM  THE  BATTERY,  BUT  GREATER  THAN  8  KM 

ALFA  FIRING  SECTION  IS  FIRING  ONE  MISSILE  AT  A  TARGET 
<40  KM  FROM  THE  BATTERY,  BUT  GREATER  THAN  8  KM 

CRUISE  MISSILE  IS  ERRONEOUSLY  IDENTIFIED  AS  A  FRIEND 
BY  IFF 

CRUISE  MISSILE  IDENTIFIED  AS  NOT  FOE  BY  SPEED  AND 
ALTITUDE  CRITERIA,  IFF  IS  NONOPERATION 

ALFA  FIRING  SECTION  IS  OPERATIONAL 

ALFA  FIRING  SECTION  IS  NONOPERATIONAL 

BATTERY  IS  NONOPERATIONAL,  NO  DETECTION  CAPABILITY 

BATTERY  IS  NONOPERATIONAL,  NO  FIRING  CAPABILITY 

ALFA  SECTION  MISSILE  "NO  KILL"  FOLLOWED  BY  RANGE  (KM) 
OF  APPROACHING  TARGET 


Table  3-3.1 
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result  printouts 


NUMBER  PRINTED  MEANING 

26  BRAVO  SECTION  MISSILE  "NO  KILL"  FOLLOWED  EV  RANGE  (KM) 

OF  APPROACHING  TARGET 

66  IFF  IS  NONOPERATIONAL 


TARGET  "KILLS"  ARE  SPELLED  OUT,  I.E.  "KILL",  FOLLOWED 
BY  THE  RANGE  FROM  THE  UNIT  AT  WHICH  THE  TARGET  WAS 
DESTROYED 

"PENETRATION"  IS  PRINTED  WHEN  A  CRUISE  MISSILE  APPROACHED 
WITHIN  S  KMs  OF  THE  BATTERY.  IN  THIS  INSTANCE  THE 
UNIT  IS  CONSIDERED  PENETRATED  AND  DESTROYED 


Table  3-3.2 
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TV.  DATA  ANALYSIS 


A.  PURPOSE 

The  purpose  of  this  chapter  is  to  demonstrate  how  the  TI  59  may  be 
used  to  analyze  data.  Rather  than  use  assumed  or  contrived  data,  an 
actual  experiment  was  conducted  by  the  authors  for  illustrative  purposes. 
Strict  requirements  for  random  sampling  were  not  met  in  conducting  the 
experiment  but,  again,  the  purpose  of  this  thesis  is  to  demonstrate  the 
capabilities  of  the  TI  59  rather  than  to  make  inferences  or  draw  broad 
conclusions  from  the  experimental  data. 

Before  discussing  methods  for  data  analysis,  the  scope,  design,  and 
methodology  of  the  experiment  will  be  presented  in  sufficient  detail  to 
make  the  data  analysis  meaningful.  Presentation  of  the  experiment  wi^l 
be  followed  by  a  detailed  discussion  of  a  TI  59  program  designed  to 
compute  measures  of  central  tendency  and  spread  for  sample  data.  Chapter  V 
discusses  a  TI  59  program  which  may  be  used  to  make  statistical  in¬ 
ferences  using  the  same  experimental  data. 

3.  FITNESS  EXPERIMENT 

1 .  Scope  of  the  Experiment 

The  experiment  was  conducted  to  test  the  hypothesis  that 
different  physical  conditioning  programs  result  in  different  levels  of 
physical  fitness.  Six  different  conditioning  programs  were  evaluated 
using  five  tests.  The  scope  of  the  experiment  was  limited  to  testing 
the  strength  and  endurance  of  selected  upper  body  muscles,  together  with 
overall  cardiovascular  fitness.  A  completely  comprehensive  fitness 
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evaluation  would  also  include  lower  body  strength  and  endurance  as  well 
as  muscular  flexibility  and  agility.  Other  factors  such  as  diet,  use  of 
alcohol,  tobacco  and  other  drugs  would  also  be  requisite.  This  experi¬ 
ment  was  limited  to  the  examination  of  thirty-six  male  subjects  by  the 
two  authors  to  determine  cardiovascular  efficiency,  bicep  strength., 
bicep  endurance,  pectoral  strength  and  pectoral  endurance.  Thirty  of 
the  subjects  were  military  officers  attending  the  Naval  Postgraduate 
School.  The  other  six  subjects  were  weightlifters  who  trained  at 
Bailey's  Gym  in  Seaside,  California.  The  subjects  varied  in  age  from 
nineteen  to  thirty-seven  but  were  predominantly  in  their  early  thirties. 
While  it  is  recognized  that  strict  requirements  for  random  sampling 
requisite  for  statistical  analysis  were  not  met  in  conducting  the  experi¬ 
ment,  it  should  be  pointed  out  that  there  is  every  reason  to  believe 
that  the  subjects  examined  were  representative  of  the  population  from 
which  students  are  continually  drawn  for  the  institution.  Strict  inference 
to  any  specified  population  will  not  be  made,  but  as  mentioned  earlier 
the  purpose  of  this  discussion  is  to  illustrate  the  use  of  programs  de¬ 
veloped  for  the  TI  59,  not  to  make  inferences  from  the  data  gathered  in 
this  particular  experiment. 

Since  the  terms,  strength,  endurance,  and  cardiovascular 
efficiency  are  subject  to  a  myriad  of  interpretations,  the  definitions 
used  for  this  analysis  are  presented  before  proceeding  to  discuss  the 
test  methodology.  Muscular  strength  is  defined  as  the  ability  to 
exert  maximum  force  against  an  object,  while  muscular  endurance  relates 
to  the  ability  to  exert  force  which  is  not  necessarily  maximal  over  an 
extended  period  of  time.  Cardiovascular  efficiency  relates  to  how 
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well  the  heart,  lungs  and  blood  vessels  work  together  without 
strain  ^Ref.  9J 

Before  analyzing  the  results  a  thorough  explanation  of  the 
experimental  design  and  methodology  is  provided. 

2.  Experimental  Design 

The  subjects  were  divided  into  six  categories  (each  consisting 
of  six  individuals)  based  on  their  exercise  programs.  The  following  six 
categories  were  used:  Category  I  -  Individuals  who  had  participated  in 
no  exercise  over  the  last  six  months;  Category  II  -  Individuals  who  ran 
more  than  thirty  miles  per  week  and  who  did  so  for  at  least  the  last 
six  months;  Category  III  -  Individuals  whose  exercise  consisted  solely 
of  lifting  weights,  whether  for  power  or  building  physique,  and  who 
lifted  at  least  twice  a  week;  Category  IV  -  Individuals  who  lifted 
weights  at  least  twice  a  week,  who  ran  twenty  or  more  miles  per  week 
and  who  did  both  the  lifting  and  running  for  at  least  the  last  six 
months;  Category  V  -  Individuals  who  lifted  weights  at  least  twice  a 
week,  who  ran  between  ten  and  nineteen  miles  per  week,  and  who  did 
the  lifting  and  the  manning  for  at  least  the  last  six  months;  and 
Category  VI  -  Individuals  who  did  not  run  or  lift  weights  but  who 
participated  in  basketball,  racquetball,  bicycling,  or  judo  on  a  regular 
basis. 

Five  tests  were  administered  to  each  subject  in  each  category. 
Cardiovascular  efficiency  was  measured  using  the  Pipes  Test  for  Cardio¬ 
vascular  Health,  which  consists  of  the  following  seven  steps: 

1.  Have  the  subject  lie  quietly  on  a  pad  for  ten  minutes. 

Take  a  pulse  reading  for  the  resting  heart  rate. 
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3.  Have  the  subject  sit  on  a  chair  with  his  arms  folded 
across  his  chest. 

4.  With  his  arms  folded,  have  the  subject  stand  up  and 
sit  down  twice  every  five  seconds  for  three  minutes. 

5.  Take  a  reading  for  the  heart  rate  immediately  after 
sitting  down  at  the  conclusion  of  the  three  minute 
exercise . 

6.  Take  readings  at  30-second,  60-second,  and  120- 
second  intervals  after  the  exercise  to  measure 
recovery  heart  rate. 

7.  Using  a  table  developed  by  Pipes  and  the  five  heart 
rate  readings  determine  the  cardiovascular  health 
score . 

This  test,  developed  by  an  exercise  physiologist,  [Itef.  9J  is 
based  on  the  premise  that  the  heart  rate  describes  an  individual's 
fitness  in  three  areas:  how  much  oxygen  he  needs,  how  much  blood  his 
heart  has  to  pump  to  meet  the  oxygen  need,  and  how  hard  the  heart  works. 
Individuals  with  a  low  level  of  fitness  do  not  extract  oxygen  efficiently 
from  the  blood  so  must  pump  more  blood,  resulting  in  a  higher  heart  rate. 
In  conducting  the  test,  the  same  examiner  took  the  pulse  readings  at 
the  different  intervals  using  the  radial  artery  in  the  subject's  wrist. 
The  heart  rate  was  monitored  for  ten  seconds  then  multiplied  by  six  to 
obtain  the  number  of  heart  beats  per  minute. 

The  muscular  strength  test  used  for  biceps  was  the  maximum 
standing  curl  the  individual  could  accomplish  using  a  barbell.  The 
proper  technique  was  demonstrated  to  each  subject  by  the  examiner  and 
lifts  where  the  individual  "cheated"  by  swinging  the  weight  or  arching 
his  back  were  not  counted.  Bicep  endurance  was  measured  by  the  number 
of  curl  repetitions  the  subject  performed  with  a  55  pound  weight.  These 
two  tests  were  predicated  on  the  generally  accepted  basis  that  curling 
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is  the  primary  bicep  exercise  and  that  low  repetition,  high  resistance 


Ref.  Jj 


Fifty-five 


exercises  are  best  for  developing  endurance 
pounds  of  weight  were  selected  prior  to  the  experiment  as  a  low 
resistance  weight  well  below  each  subject's  strength  capability  and 
therefore  in  compliance  with,  expert  opinion  that  the  force  used  for 
endurance  testing  should  be  considerably  below  the  individual's  static 
force  capability  [lief.  1J 


In  a  similar  fashion,  the  bench  press  was  used  to  test  pectoral 
strength  and  endurance.  Pectoral  strength  was  measured  by  the  maximum 
weight  that  the  individual  could  bench  press,  while  pectoral  endurance 
was  measured  by  the  number  of  bench  press  repetitions  he  could  perform. 

Each  subject  was  tested  using  an  identical  sequence  of  events. 
Initially  each  person  was  given  an  instruction  sheet  (table  B-l  in 
Appendix  B)  which  explained  the  purpose  of  the  experiment  and  defined 
those  attributes  to  be  measured,  i.e.  muscular  endurance,  muscular 
strength,  and  cardiovascular  health.  The  subject  was  then  asked  to 
complete  a  questionnaire  (table  3-2,  Appendix  B)  concerning  certain 
aspects  of  his  medical  history.  Each  man  was  instructed  to  stop  the 
testing  if  he  felt  any  significant  level  of  pain.  He  was  then  asked 
to  complete  a  form  disclosing  his  name,  age,  weight,  and  height 
(table  B-3,  Appendix  B) .  Next,  one  of  the  examiners  questioned  the 
subject  concerning  his  exercise  program  over  the  last  six  months  and 
made  a  subjective  judgement  as  to  which  of  the  six  categories  he 
belonged  in. 

After  this  administrative  procedure  was  completed,  the  actual 


testing  was  begun  with  the  bench  press  test  described  above.  (The 
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same  examiner  tested  each  of  the  thirty-six  subjects  in  both  the  strength 
and  endurance  exercises  in  order  to  minimize  any  variance  due  to  test 
administration.)  The  subject  was  shown  how  to  do  the  bench  press, 
allowed  to  practice  once  if  desired,  and  then  tested  for  the  maximum 
number  of  repetitions  he  could  perform  with  100  pounds.  The  number 
of  repetitions  was  recorded  and  the  subject  was  allowed  a  three  minute 
rest  before  being  tested  for  his  maximum  bench  press.  The  examiner 
estimated  the  amount  that  each  subject  could  bench  press  and  set  up 
the  weights  accordingly.  All  adjusting  of  weight  was  done  by  the 
examiner  so  that  the  subject's  lift  capability  was  not  degraded.  The 
man  was  then  asked  to  bench  press  the  weight  set  up  for  him.  If  he 
was  able  to  make  the  lift  ten  pounds  were  added  and  he  was  asked  to 
try  again.  If  he  failed  the  second  attempt  he  was  given  credit  for 
five  pounds  less  than  he  attempted.  For  example,  if  an  initial  attempt 
of  165  pounds  was  successful  and  a  subsequent  attempt  of  175  pounds  was 
missed,  then  the  score  was  recorded  as  170  pounds. 

After  another  three  minute  rest  the  subject  was  tested  on  the 
number  of  times  he  could  curl  55  pounds.  This  was  followed  by  another 
three  minute  rest  before  testing  for  his  maxi  mum  curl  capability.  Once 
again,  the  man's  maximum  lift  was  estimated  by  the  examiner  and  all 
adjusting  of  the  weight  was  done  by  the  examiner. 

Following  the  four  lift  tests,  the  Pipes  Cardiovascular  Health 
Test  [nef.  9^  was  administered  by  the  second  examiner  and  the  subject's 
testing  was  completed.  Before  examining  the  test  results,  a  procedure 
for  scoring  the  tests  was  requisite.  Accordingly,  the  scoring  procedure 
explained  in  the  following  discussion  was  decided  upon. 
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3.  Scoring  Methodology 


The  heart  rates  recorded  during  the  Pipes  Cardiovascular  Test 

were  scored  using  table  B-4  in  Appendix  B.  Each  subject's  score  for 

resting  heart  rate,  heart  rate  immediately  after  the  exercise  and  heart 
rate  at  the  30-second,  60-second,  and  120-second  intervals  was  aggregated 
to  a  total  score  ranging  from  zero  to  one  hundred.  This  score  was  then 
used  as  the  measure  of  cardiovascular  fitness  for  comparative  analysis. 

In  order  to  compare  muscle  strength  among  the  subjects  and 
among  the  categories,  it  was  necessary  to  adjust  each  subject's  lift  for 
varying  sizes  and  body  structures.  Accordingly,  each  roan's  maximum  curl 
and  maximum  bench  press  were  divided  by  his  body  weight,  resulting  in 

an  adjusted  score  for  each  lift.  These  two  adjusted  scores  were  then 

added  together  to  yield  an  upper  body  strength  measure.  For  example, 
let,  Sj_3  be  the  strength  measure  for  the  third  subject  in  Category  I 
(where  the  first  subscript  indicates  the  category  and  the  second  indicates 
the  subject  within  the  category) .  The  following  formula  may  then  be  used 
to  obtain  the  strength  score  for  the  third  subject  tested  in  Category  I: 
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maximum  bench  +  maximum  curl 
body  weight 


As  discussed  earlier,  muscle  endurance  was  measured  for  the 
same  two  areas  tested  for  strength  -  the  biceps  and  the  pectorals.  The 
bicep  endurance  was  measured  by  the  maximum  number  of  curl  repetitions 
performed  with  55  pounds,  while  the  pectoral  endurance  was  measured  by 
the  maximum  number  of  bench  presses  accomplished  with  100  pounds.  As 


in  the  case  of  strength,  an  adjustment  was  made  for  the  subject's  body 
weight.  In  the  case  of  endurance,  however,  the  amount  of  weight  lifted 
(which  was  55  pounds  for  the  curl  and  100  pounds  for  the  bench  pressl 
was  divided  by  the  subject's  body  weight  and  then  multiplied  by  the 
corresponding  number  of  repetitions  lifted.  These  two  scores  were  then 
summed  as  the  endurance  index.  For  example,  let  be  the  endurance 
score  for  the  third  subject  tested  in  Category  I.  The  following  formula 
then  obtains: 
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_ 100  lbs _ 

body  weight  of  subject 


♦(number  of  bench  press 
repetitions! 


55  lbs 


body  weight  of  subject 


(number  of  curl 
repetitions) 


4.  Test  Results 

Tables  B-5  through  B-10  in  Appendix  B  reflect  the  results  of  the 
experiment  for  each  of  the  six  categories  tested.  For  example,  table  B-5 
depicts  the  age,  weight,  cardiovascular  scare,  adjusted  strength  score 
and  adjusted  endurance  score  for  each  of  the  six  subjects  tested  in 
Category  I.  Appropriate  references  are  made  at  table  B-5  for  the  develop¬ 
ment  of  the  final  cardiovascular,  strength  and  endurance  scores.  The 
scoring  methodology  section  of  the  chapter  provides  a  detailed  explanation 
of  the  rationale  and  methodology  for  deriving  these  scores. 

C.  TX  59  FROGRAM  FOR  DATA  ANALYSIS 


Having  developed  the  experimental  design,  the  scoring  methodology, 
and  the  test  results  it  is  now  possible  to  analyze  the  data-  Measures  bf 


central  tendency  and  spread  will  be  used  to  illustrate  an  application  of 
the  TI  59  in  analyzing  sample  data.  The  measures  of  central  tendency 
used  for  this  illustration  are  the  mean,  "x  ,  and  the  median.  The 
measures  of  spread  used  are:  the  sample  variance,  s“;  the  standard 
deviation,  s  ;  the  mean  absolute  deviation,  MAD:  the  mean  squared 
deviation,  MSD;  the  root  mean  squared  deviation,  RMS D ;  and  the  range. 

A  TI  59  program  will  now  be  described  in  detail  which  computes  these 
measures  followed  by  an  example  applying  the  program  to  the  results  of 
the  fitness  experiment. 

1 .  TI  59  Capabilities 

The  TI  59  has  been  hard-wired  to  calculate  the  sample  mean  and 
variance  as  well  as  MSD.  As  described  in  the  TI  59  Personal  Programming 
Manual  [Ref .  $  if  each  datum  is  entered  into  the  calculator  followed 
by  pressing  the  £+  key,  the  calculator  will  sum  each  data  entry,  :cj_, 
into  register  one,  sum  the  squares  of  xj_  into  register  two,  store  the 
number  of  data  entries  in  register  three,  and  calculate  x,  s2,  and  MSD. 
(By  definition,  s2  =  1/n-l  £(x^  -  x) 2  is  the  unbiased  estimator  for 
T 1  while  MSD  =  (n-l)s2/n  is  the  maximum  likelihood  estimator.) 

Pressing  the  key  will  yield  the  mean,  INV  "x  will  display  s2  and  2nd  Op 
11  will  display  MSD.  If  these  are  the  only  measures  desired  then 
utilization  of  the  £+•  key  is  the  most  expedient  method  of  obtaining 
them.  The  TI  59  statistics  module  has  a  program  (Program  03)  which 
computes  these  same  measures  as  well  as  the  middle  value  (MIDVAL)  . 
Additionally,  Program  03  stores  each  data  entry  beginning  with  register 
31.  Program  03  also  computes  a  number  of  other  quantities  not  germane 


to  an  analysis  of  the  data  gathered  in  the  experiment  discussed  pre¬ 


viously  in  this  chapter.  Since  this  results  in  a  slightly  longer  run 
time  for  each  computation,  a  program  has  been  written  by  the  authors 
which  exploits  the  hard-wire  capabilities  of  the  £+  key,  computes  MAD 
and  range  in  addition  to  the  other  measures  discussed,  and  stores  the 
data  for  recall  or  transformation  if  desired.  In  addition,  this  pro¬ 
gram  may  be  used  with  the  TI  59  Master  Module  if  the  Statistics  Module 
is  not  available.  The  following  section  describes  the  program  in  detail. 

1 .  Univariate  Data  Program 

In  order  to  facilitate  the  description  of  this  program  a  flow¬ 
chart  (figure  4-1),  has  been  included  at  the  end  of  the  chapter. 

Comments  in  the  paper  are  keyed  to  figure  4-1  by  numbered  circles  for 
easy  reference. 

The  program  is  initialized  by  pressing  E'  (figure  4-1, (T)). 
Initialization  entails  clearing  all  of  the  data  registers,  lowering 
flag  0  the  purpose  of  which  will  be  addressed  later,  and  storing  31  in 
register  iO .  Register  30  is  used  as  a  post  office  for  indirect  address¬ 
ing.  In  this  particular  program  this  means  that  data  are  stored  in  the 
register  indicated  by  register  30.  For  example,  after  initialization, 
register  30  contains  31.  The  sequence  x-^,  STO  2nd  IND  30  will  result  in 
x^  being  stored  in  register  31.  When  the  initialization  routine,  2nd  E', 
is  complete  the  display  will  contain  the  value  31.  Each  datum  may  now  be 
entered  successively  followed  by  pressmc  A.  The  routine  at  Label  A 
begins  by  storing  x^  in  a  working,  register  13  (figure  4-1,^^)  . 

Xj_  is  then  stored  permanently  beginning  in  register  31.  is  stored 


in  register  31,  then  register  30  the  indirect  storage  address  is  incre¬ 


mented  by  1  so  that  x-,  will  be  stored  in  register  32,  X3  in  register  3  3 

and  xn  in  register  31  +  n  -  1.  A  total  of  sixty-nine  entries  mav  be 

made  using  registers  31  through  99  for  data  storage.  Registers  0 

through  29  are  used  to  make  the  requisite  computations  of  central 

tendency  and  spread.  After  each  datum,  x^,  has  been  stored,  the  program 

checks  to  see  if  flag  0  is  raised  (figure  4-1, ^7))  .  If  flag  0  is 

raised  this  indicates  that  a  data  entry  has  been  made  previously,  ie . , 

the  current  x^  is  not  x, .  In  this  event  the  program  skips  to  Label  x. 

If  flag  0  is  not  raised,  ie.,  the  current  x^  is  x^,  then  is  recalled 

from  the  working  register,  register  18,  and  stored  in  register  12  as 

the  minimum  x^  and  register  13  as  the  maximum  x^.  Future  entries  may 

then  be  checked  against  register  12  to  determine  which  value  is  lower. 

If  a  current  Xj_  is  lower  than  the  value  in  register  12  then  it  will 

replace  it  as  x^,.^.  Similarly,  subsequent  entries  may  be  checked 

against  register  13  in  order  to  retain  x-,^-  After  storing  x-  in 

register  12  and  in  register  13  the  program  internallv  calls  the  key 

which,  as  disscussed  previously,  will  sum  x^  into  register  01,  sum 

x,“  into  register  02,  sum  the  number  cf  entries  into  register  03  and 

compute  x”,  s'-,  and  MSO .  Flag  0  is  then  raised  so  that  subsequent 

entries  will  skip  to  Label  x  and  reciace  x_,„  cr  x  as  aoorooriate . 

-  mm  max 

The  program  then  recalls  the  number  of  entiles,  n,  into  the  display  and 
stops  awaiting  the  next  entry  'figure  4-1,  (7^)  .  The  second  entry  x-,  will 
now  be  stored  temporarily  in  the  working  register,  register  18,  and 
permanently  in  register  32.  The  indirect  addressina  register,  register 
30,  is  incremented  by  1  for  the  next  entry  and  the  program  then  checks 


to  see  if  flag  0  is  raised.  Since  this  is  not  the  first  entry,  the  flag 
will  be  raised  causing  the  program  to  skip  to  Label  x  (figure  4-1,  ^5^)  ). 
The  first  step  under  Label  x  is  to  recall  from  register  12  and 

store  it  in  the  test  register  R^.  For  this  particular  iteration,  x^ 
will  be  in  register  12  since  the  first  entry  was  both  the  maximum  and 
the  minimum  value  processed  as  described  earlier.  The  program  then 
recalls  the  current  x^  (x2  in  this  instance)  from  the  working  register, 
register  18.  The  display  value,  x, ,  is  checked  against  the  R™  value, 
x^,  to  see  if  the  display  value  is  less  than  the  value.  If  so  the 
program  skips  to  step  57  where  x-,  is  stored  in  register  12,  replacing 
Xt_  as  the  lowest  data  entry  (figure  4-1, (*T))  .  If  the  display  value  is 
not  less  than  the  value  then  the  program  recalls  Xj,^  from  and 
stores  it  in  Rm.  x^  is  recalled  from  the  working  register,  register  13, 
into  the  display.  This  time  the  program  checks  to  see  if  the  display 
value,  x^  is  greater  than  the  R^  value,  x^^.  If  so,  the  program  skips 
to  step  62  where  x^  is  stored  in  register  13  as  the  new  x^^  (fig.  4-1, 
©>  .  The  program  then  computes  the  MIDVAL  by  recalling  from 

register  12  and  ;tnax  from  register  13,  summing  them  and  dividing  by 
2,  and  storing  in  register  14.  Next  the  range  is  computed  by  sub¬ 
tracting  x^^  from  x^a^.  The  range  value  is  stored  in  register  15 
(figure  4-1,^)].  The  program  then  loops  back  to  the  +  key  to 

compute  the  mean,  variance  and  MSD,  (figure  4-1,^^).  The  number  of 
entries  is  recalled  and  displayed  awaiting  the  next  entry.  This  process 
is  repeated  until  all  of  the  data  have  been  processed. 

The  outputs  of  the  program  may  be  recalled  as  shown  in  table  3-11 
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in  Appendix  3.  The  mean  is  displayed  by  pressing  x",  the  variance  by 
pressing  Inv  x",  and  MSD  by  pressing  2nd  op  11.  The  lowest  data  point, 
xmin '  may  136  discovered  by  recalling  register  12  while  the  highest  data 
entry,  x^^,  may  be  recalled  from  register  13.  Recalling  register  14 
will  display  the  MIDVAL  and  the  range  may  be  found  by  recalling 
register  15.  Each  of  the  original  data  entries  may  be  recalled  if 
desired  beginning  with  x^_  in  register  31.  MAD  is  computed  by  pressing 
2nd  A1  which  calls  a  different  subroutine.  This  subroutine  recalls  x 
which  was  computed  under  Label  A  and  stores  it  in  Register  16.  The 
number  of  entries,  n,  is  recalled  from  register  3  and  stored  in 
register  7  to  be  used  as  a  decrement  register.  Register  20  contains 
31  which  is  used  to  indirectly  address  the  datum  which  have  been  stored 
beginning  with  register  31.  The  program  recalls  each  x^  using  register 
2C  and  subtracts  x\  The  absolute  value  of  the  difference  is  summed 
into  register  19.  The  program  does  this  successively  for  each  x^  until 
the  decrement  register,  register  20,  is  equal  to  zero  indicating  that 
each  x^  has  been  processed.  The  sum  of  the  absolute  values  of  the 
deviations  from  the  mean  is  recalled  from  register  19  and  divided  by 
n  which  is  recalled  from  register  3.  This  value,  the  mean  absolute 
deviation  is  displayed  completing  the  subroutine  A'  processing.  All 
the  values  discussed  earlier  are  still  intact  and  may  be  recalled  if 
needed.  Table  3-12,  Appendix  B,  is  a  program  listing  for  the  uni¬ 
variate  program. 

0.  APPLICATION  OF  THE  TI  59  UNIVARIATE  PROGRAM 

The  cardiovascular  scores  for  Category  I  provide  a  ready  example  for 
the  use  of  the  univariate  program  to  calculate  measures  of  central 


tendency  and  spread.  After  the  program  card  has  been  read  in,  the 


program  is  initialized  by  pressing  2nd  S' .  The  cardiovascular  scores 
for  Category  I  (table  3-5)  are  entered  into  the  calculator  as  follows: 
56.5,  A;  58,  A;  58,  A;  44.5,  A;  33.5,  A;  40,  A.  The  instructions 
contained  in  table  B-li  may  then  be  used  to  obtain  the  desired  statistics 
For  this  particular  example:  2nd  x  yields  the  mean,  48.4;  RCL  14  dis¬ 
plays  the  MIDV.AL,  45.75;  RCL  15  displays  the  range,  24.5;  2nd  A*  yields 
the  mean  absolute  deviation,  9;  2nd  Op  11  displays  the  mean  squared 
deviation,  92.9;  2nd  Op  11  ,~\l  x  ,  calculates  the  root  mean  squared 

deviation,  9.6;  INV  2nd  x  recalls  the  standard  deviation,  10.56;  INV 

_  -> 

2nd  x,  x^  calculates  the  variance  111.5;  RCL  12  displays  the  lowest  data 
entry,  33.5  and  RCL  13  displays  the  highest  data  entry,  58.  To 
calculate  the  sample  statistics  for  another  category  or  for  a  different 
test  the  user  need  only  push  2nd  S’  to  re-initialize  and  then  enter  the 
relevant  data.  Statistics  have  been  calculated  for  the  age,  weight, 
cardiovascular  scores,  endurance  scores,  and  strength  scores  for  each 
of  the  six  categories.  Tables  4-1  through  4-6  display  these  statistics. 
Rather  than  discuss  each  of  these  tables  in  depth,  one  example  is  pro¬ 
vided  relative  to  the  interpretation  of  the  sample  statistics. 

The  cardiovascular  mean  for  Category  I,  48.4,  indicates  average 
cardiovascular  fitness  using  the  Pipe's  test  which  is  based  on  a  scale 
from  0  to  100.  Three  measures  of  spread i (standard  deviation  10.56, 
root  mean  squared  deviation,  9.6,  and  mean  absolute  deviation,  9)  are 
approximately  equal  to  ten,  a  rather  high  variability  in  this  case. 

The  range,  24.5,  also  indicates  that  the  data  are  quite  spread  out. 
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x  .  of  33.5  and  x„  of  58,  the  bounds  of  the  sample  data  indicate  that 
min  max  T 

the  cardiovascular  fitness  of  sedentary  people  varies  from  poor  to 
average . 

Inferences,  subject  to  the  sampling  limitations  already  discussed, 
may  also  be  made  about  the  strength  or  endurance  of  sedentary  people 
using  the  data  from  table  B-5 .  Similarly,  the  statistics  for  the 
other  categories  may  be  used  to  make  inferences  about  the  strength, 
endurance,  or  cardiovascular  fitness  of  those  who  run  over  thirty  miles 
per  week  (Category  II)  or  those  who  lift  weights  (Category  III)  or 
any  cf  the  other  three  categories.  Programs  have  also  been  written  for 
the  TI  59  which  allow  a  user  to  develop  confidence  intervals  for  these 
sample  statistics  jRef.  14^  . 

The  next  chapter  will  discuss  a  program  for  one  factor  analysis 
of  variance  and  then  apply  the  program  to  the  fitness  data  to  illustrate 
statistical  inference  with  the  TI  59. 


71 


Figure  4-1.4 
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V.  STATISTICAL  INFERENCE 

A.  PURPOSE 

The  purpose  of  this  chapter  is  to  illustrate  a  method  of  statistical 
inference  using  the  TI  59.  As  in  chapter  IV  the  intent  is  to  demonstrate 
the  capabilities  of  the  TI  59  rather  than  to  emphasize  statistical  prin¬ 
ciples.  One-way  analysis  of  variance (ANOVA)  will  be  used  to  demonstrate 
statistical  inference  using  the  fitness  data  discusser  in  the  preceeding 
chapter.  The  variations  of  the  underlying  populations  represented  by  the 
six  categories  are  assumed  to  be  unknown  but  equal  for  this  illustration. 

B.  TI  59  PROGRAM  FOR  ANALYSIS  OF  VARIANCE 

In  testing  the  hypothesis  that  the  population  means  for  each  of  the 
six  test  categories  are  equal,  HQ  is  typically  rejected  if  the  F  ratio 
exceeds  the  critical  F  value  in  a  standard  table  for  the  desired  test 
level  (typically  5%).  Alternatively,  using  the  TI  59  to  its  full 
advantage,  prob-value  may  be  used  to  test  Hg .  Prob-value  is  a  method  of 
testing  whether  or  not  the  null  hypothesis  is  supported  by  the  data.  In 
the  case  of  the  F  ratio,  prob-value  is  the  probability  that  the  F  ratio 
would  be  as  large  or  larger  than  the  value  actually  observed  if  Hq  were 
true.  This  is  the  right  hand  tail  area,  Q  (f) ,  where 
Q(f)  =  Pr  (F  >  f) 

Prob-value  has  the  advantage  that  analysis  is  not  restricted  to  arbi¬ 
trarily  established  test  levels  such  as  5%  or  10%  or  to  use  of  standard 
published  tables.  The  TI  59  Statistics  Module  has  an  F  distribution 
program  (Program  22)  which  computes  the  tail  area  of  an  F  curve  where 


the  curve  is  defined  by  the  degrees  of  freedom  in  the  numerator  and  the 


denominator.  A  series  expansion  is  used  to  approximate  the  integral  to 
determine  Q(f)  [Ref.  lcj  .  If  HQ  is  true,  indicating  that  all  of  the 
observations  are  from  the  same  normal  population,  then  the  prob-value, 
2(f),  will  be  large.  Conversely,  if  is  false  then  the  prob-value 
will  be  small.  If  the  prob-value  is  sufficiently  small  (as  determined 
by  the  decision-maker)  then  is  rejected  and  the  conclusion  is  formed 
that  there  must  be  a  difference  ir.  the  population  means  somewhere 


While  a  classical  test  or  prob-value  may  facilitate  rejection  or 
acceptance  of  the  null  hypothesis,  no  insight  is  provided  as  to  which 
means  differ,  given  that  is  rejected.  There  is  an  efficient  method 
developed  by  Sheffe  §ef.  l]  for  computing  confidence  intervals  for 
the  difference  between  means.  If  the  physical  fitness  example  discussed 
earlier  is  used,  then  Sheffe 's  development  may  be  used  to  make  the 
following  statements  with  95%  confidence: 


(u,  -  ut)  -  (x  -  Xj  ±  -V (r-i)£t  Sp-V-A,  +  7T, 

±  _V(r-i)F„s^/  n, ~  n3 
=  (*-*.)  ±  ^r-O^p^n/  7T3 

Where  F  is  the  critical  value  of  F  which  leaves  5%  in  the  upper  tail, 


Sp  is  the  square  root  of  the  pooled  variance,  r  is  the  number  of  means 


compared,  and  n  is  the  size  of  each  of  the  samples.  In  the  fitness 
example  this  equates  to  making  confidence  statements  about  the  difference 
in  fitness  between  the  six  categories  where  r  is  six  and  n  is  6  for  each 
of  the  samples.  To  facilitate  multiple  comparisons  a  contrast  of  means 
is  used.  This  contrast  may  be  written  as: 

E  CA/.  where  [c,  =  0 

It  is  then  possible  to  develop  the  following  formula  which  includes  all 
possible  contrasts  with  95%  confidence: 

EciM=  Eqx;  ±-x/F)?s  SpVE^f 

If  the  value  0  is  included  in  a  confidence  interval  then  there  is  no 
basis  for  believing  that  the  population  means  differ  while  if  0  is  not 
included  then  the  conclusion  is  drawn  that  the  means  do  in  fact  differ. 

A  program  has  been  developed  by  Dr.  P.W.  Zehna  jnef.  l^j  for  the 
TI  59  which  computes  the  elements  of  an  ANOVA  table  (table  B-15)to 
include  the  F  ratio  and  prob-value  discussed  above.  Basically  the 
program  exploits  the  TI  59  F  distribution  program  for  determining 
prob-value  (Program  22)  after  using  Program  15  of  the  Statistics  Module 
[Ref.  loj  to  calculate  the  F  statistic.  The  program  then  uses  Scheffe's 
multiple  contrasts  to  determine  which  population  means  differ  given 
that  the  null  hypothesis  is  rejected. 

A  flowchart  (figure  5-1) ,  user  instructions  (table  B-13) ,  and  a  listing 
of  the  actual  program  steps  (table  B-14)  are  provided  to  facilitate 
description  of  the  ANOVA  program.  The  program  takes  data  input  by  rows 
and  outputs  the  elements  of  an  ANOVA  table  (table  B-15)  sequentially 


as  indicated  by  the  number  in  each  block  of  the  table.  A  row  of  data 
constitutes  a  sample  as  in  the  example  at  table  B-16 . 

The  program  begins  by  using  Program  06  of  the  Statistics  Module  to 

enter  the  data.  After  initialization  with  2nd  E'  each  x. .  is  entered 

JO 

followed  by  pressing  Label  A.  When  one  complete  row  has  been  entered  a 
press  of  2nd  B'  causes  the  calculator  to  compute  the  row  or  sample  x  and 
a  press  of  2nd  C'  results  in  computation  of  the  MSD .  These  two  steps 
must  be  performed  after  each  row  has  been  input  so  that  the  calculator 
will  know  when  a  new  row  is  being  entered.  When  all  of  the  data  have 
been  entered  using  this  scheme  (table  3-13)  the  sequence  RST,  A  begins 
the  ANOVA  Table  calculations.  The  first  step  under  Label  A  is  to  call 
Program  15  of  the  TI  59  Statistics  module  which  computes  the  F  ratio. 

In  the  process  of  computing  the  F  ratio  the  other  elements  of  the  ANOVA 
table  (table  B-15)  are  computed  and  stored  except  for  the  prob-value . 

To  fill  in  the  values  for  the  ANOVA  table  all  that  is  required  is 
successive  pushes  of  R/S  as  indicated  in  table  B-13.  For  example,  the 
first  R/S  displays  the  degrees  of  freedom  for  the  numerator  while  the 
fourth  R/S  displays  the  degrees  of  freedom  for  the  denominator.  The 
program  essentially  recalls  and  displays  the  calculations  of  Program  15 
of  the  Statistics  Module  to  build  the  ANOVA  table .  To  compute  prob- 
value  the  program  internally  calls  Program  22  of  the  Statistics  Module. 
The  user  need  only  press  R/S  as  indicated  in  table  B-13  which  causes  the 
calculator  to  recall  the  degrees  of  freedom  for  the  numerator  and  the 
denominator  used  in  Pgm  15  and  transfer  them  to  Program  22  to  define 
the  F  Distribution.  The  F  statistic  calculated  in  Program  15  is  then 
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recalled  and  transfered  to  Program  22  resulting  in  o(f),  the  probability 
that  F  >  f.  This  prob-value  may  then  be  used  to  accept  or  reject  the 
null  hypothesis.  If  the  null  hypothesis  is  accepted  then  the  analysis 
is  completed.  However,  if  HQ  is  rejected,  the  next  step  entail?  the  use 
of  Scheffe's  contrasts  to  determine  which  means  differ. 

To  use  the  ANOVA  program  (table  B-14)  for  posterior  contrasts  with 
Scheffe’s  formulas  the  user  initializes  the  routine  by  pressing  2nd  E' . 

Then  c^ ,  x^  and  are  entered  for  each  row  as  shown  in  table  3-13.  The 

Cj_ '  s  are  the  coefficients  used  to  determine  which  means  are  contrasted 
as  discussed  previously.  To  contrast  u-j_  and  U2/  Cj.  =  1,  c,  =  -1  and  all 
other  c^  =  0.  To  contrast  u->  and  u^ ,  C2  =  1,  C3  =  -1  and  all  other 
c^  =  0 .  As  these  data  are  input,  the  program  uses  a  'loop'  to  calculate 
which  is  stored  in  register  03  and  £  Cj  X]  which  is  stored 
in  register  06,  (figure  5-1).  Register  04  is  used  as  a  counter  to 
display  the  number  of  row  entries.  After  each  c^_ ,  x^  and  n^  entry, 
register  04  is  incremented  by  one.  The  program  then  transfers  to 
Label  x“ ,  displays  the  running  count  of  row  entries  and  stops  pending  the 

next  entry.  After  every  c^,  x^  and  n^  have  been  processed,  a  critical 

value  of  F  with  degrees  of  freedom  r-1  and  n(r-l)  is  entered  for  the 
desired  test  level  followed  by  2nd  A' .  The  program  recalls  the  degrees 


of  freedom  for  the  numerator,  r-1,  from  register  14  and  the  pooled 
2 

variance  Sn  from  register  29.  The  product  (r-1)  F->)S  is  formed  and 
P  _  P 


multiplied  by  the  contents  of  register  03 


Z(r»)  • 


The  square  root 


of  this  product  is  stored  in  register  05.  This  value  is  then  added  to  and 


subtracted  from  the  contents  of  register  06,  Xcixi'  to  f0Em  t^ie  desired 
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confidence  interval.  The  lower  bound  is  displayed  after  the  use  of 
2nd  A1  and  the  upper  bound  may  be  recalled  by  pressing  R/S  (table  3-13). 

An  example  wil_  now  be  provided  using  this  program  to  test  for 
differences  between  population  means  for  the  fitness  experiments. 

C.  APPLICATION  OF  THE  TI  59  ANOVA  PROGRAM 

The  null  hypothesis,  H  ,  may  be  stated  as  -  there  is  no  difference 
in  the  cardiovascular  fitness  of  those  who  do  no  exercise  (Category  I) , 
those  who  run  in  excess  of  thirty  miles  per  week  (Category  II) ,  those 
whose  exercise  consists  solely  of  lifting  weights  (Category  III)  ,  those 
who  lift  weights  and  run  in  excess  of  twenty  miles  per  week  (Category  IV) , 
those  who  lift  weights  and  run  between  ten  and  nineteen  miles  per  week 
(Category  V) ,  and  those  who  do  not  run  or  lift  weights  but  participate  in 
other  activities  such  as  basketball,  racquetball  or  bicycling  (Category 
VI) .  Table  B-35  reflects  the  cardiovascular  score  for  each  of  the  thirty- 
six  subjects  tested  by  category  as  well  as  the  mean  for  each  category. 

The  null  hypothesis  that  there  is  no  difference  between  these  category 
means  may  be  tested  using  the  ANOVA  program  with  each  of  the  categories 
constituting  a  row  for  input. 

After  the  program  card  has  been  read  in,  the  ANOVA  program  is  used 
by  first  calling  program  06  of  the  statistics  module  to  enter  the 
data.  After  initialization  with  2nd  E' ,  the  data  for  each  row  are 
entered  followed  by  A.  For  the  cardiovascular  scores  (table  3-35)  the 
first  row.  Category  I,  would  be  entered  as  follows:  56.5,  A;  58,  A; 

44.5,  A;  33.5,  A;  40,  A.  Once  the  row  data  have  been  entered  2nd  B’ 
is  pressed  to  display  the  row  mean,  48.4,  followed  by  C'  which  displays 
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the  row  MSD,  92.9.  (The  row  mean  must  be  recorded  for  use  in  posterior 
contrasts) .  The  data  are  then  entered  in  a  similar  fashion  for  the 
remaining  five  rows  (Categories  XI  through  VI) .  Once  all  of  the  data 
have  been  entered,  RST  is  pressed  to  return  the  calculator  pointer  to 
the  ANOVA  program.  A  is  then  pressed  resulting  in  calculation  of  the 
ANOT,A  Table  entries.  The  ANOVA  entries  are  recalled  with  sequential 
presses  of  R/S.  Table  B-13  discussed  earlier  contains  detailed 
instructions  on  the  use  of  the  ANOVA  program.  Table  5-1  depicts  the 
ANOVA  calculations  for  the  cardiovascular  scores  of  the  six  fitness 
test  categories.  The  prob-value  of  .00027  is  sufficiently  small  to 
cast  doubt  upon  the  null  hypothesis  that  there  is  no  difference  in  the 
cardiovascular  fitness  among  the  six  categories  tested. 

As  discussed  previously,  the  prob-value  tells  how  credible  the  null 
hypothesis  is  but  it  does  not  tell  which  categories  differ  given  that 
there  is  cause  to  reject  Hg.  However,  confidence  intervals  may  be 
established  for  contrasts  between  the  categories  using  Scheffe's 
formula  which  is  programmed  in  the  ANOVA  routine  |Ref.  13j  .  For  the 

cardiovascular  example,  the  cardiovascular  fitness  of  the  sedentary 
subjects  (Category  I)  may  be  contrasted  to  the  cardiovascular  fitness 
of  the  runner  (Category  II)  as  a  demonstration  of  the  program.  2nd  E' 
is  pressed  to  initialize  the  contrast  routine  followed  by  Cj_,  x^  and 
n^  for  each  of  the  two  rows.  For  Categories  I  and  II  the  entries  are: 


1 

R/S 

48.4 

R/S 

6 

R/S 

1 

R/S 

83.4 

R/S 

6 

R/S 
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The  appropriate  F  percentile  is  entered  followed  by  A'  to  generate 
the  desired  confidence  interval.  To  display  a  95%  confidence  interval 
for  the  difference  in  cardiovascular  fitness  between  Categories  I  and 
XI  an  F  percentile  of  2.53  (where  there  are  five  degrees  of  freedom  in 
the  numerator  and  thirty  degrees  of  freedom  in  the  denominator)  is  used 
resulting  in  an  interval  from  -65.5  to  -41.4.  Since  0  is  not  included 
in  the  interval  it  is  reasonable  to  conclude  that  there  is  a  difference 
in  the  cardiovascular  fitness  of  the  two  categories.  Table  5-2  con¬ 
tains  the  results  of  contrasting  each  of  the  six  fitness  categories. 
Only  four  contrasts  result  in  the  conclusion  that  there  is  a  difference 
between  the  categories  with  95%  confidence:  Category  I  -  Category  II 
(-65.5,  -41.4);  Category  I  -  Category  IV  (-65.5,  -4.4);  Category  II  - 
Category  III  (.6,  61.8)  and  Category  III  -  Category  IV  (-61.7,  -6.Q). 
These  results  indicate  with  95%  confidence  that  there  is  a  difference 
in  the  cardiovascular  fitness  of  those  who  run  more  than  twenty  miles 
per  week  (Categories  II  and  IV)  and  those  who  do  no  running  at  all 
(Categories  I  and  III) ,  at  least  for  those  subjects  examined. 

The  ANOVA  program  has  also  been  applied  to  the  strength  and 
cardiovascular  scores  resulting  from  the  experiment.  Table  B-36, 
Appendix  3,  reflects  the  strength  scores  of  each  of  the  thirty-six 
subjects  by  category.  The  ANOVA  results  are  contained  in  Table  5-3. 

The  prob-value  of  .12  x  10  indicates  that  the  null  hypothesis  that 
there  is  no  difference  in  the  strength  of  members  of  the  different 
categories  should  be  rejected.  Further  analysis  with  posterior  con¬ 
trasts  is  necessary  to  see  which  categories  differ.  Table  5-4  contains 
the  results  of  posterior  contrasts  with  an  F  percentile  of  2.53  for  95% 


p 


confidence  with  five  and  thirty  degrees  of  freedom  for  the  numerator  and 
denominator  respectively.  Unlike  the  cardiovascular  contrasts  there 
are  a  number  of  differences  in  the  strength  results.  Categories  III 
and  IV,  which  were  composed  of  the  most  ardent  weightlif ters ,  differs 
from  Categories  I  and  III  but  not  VI.  These  results  are  not  surprising 
in  that  they  confirm  the  hypothesis,  that  different  training  programs 
result  in  different  levels  of  fitness.  In  this  instance  where  fitness 
is  defined  as  strength,  those  who  trained  for  strength  were  in  fact 
stronger  than  those  who  did  not.  Again,  without  attempting  inference 
to  a  larger  population,  these  results  may  be  used  to  gain  insight  into 
the  probable  differences  that  might  be  tested  in  a  more  appropriately 
designed  experiment. 

Table  B-37,  Appendix  B  depicts  the  endurance  scores  for  the  thirty- 
six  subjects  by  category.  Table  5-5  reflects  the  results  of  using  the 
ANOVA  program  with  the  endurance  scores  as  input.  Once  again,  the 
prob-value  of  .137  x  10-6  indicates  rejection  of  the  null  hypothesis 
(in  this  case  that  there  is  no  difference  in  the  upper  body  endurance 
of  the  members  of  the  six  different  categories) .  The  posterior  con¬ 
trasts  (table  5-6)  indicate  that  the  weightlifters  (Categories  III, 

IV,  V)  differ  from  the  non-weightlif ters  (Categories  I,  II,  VI)  in 
upper  body  endurance  with  95%  confidence.  This  also  supports  the 
hypothesis  that  different  training  programs  result  in  different  levels 
of  fitness,  subject  again  to  the  sampling  restrictions  previously 
discussed. 
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D.  SUMMARY 


While  the  results  of  the  fitness  experiment  are  interesting,  the 
purpose  of  this  analysis  has  been  to  demonstrate  a  statistical  applica¬ 
tion  of  the  TI  59  and  not  draw  inference  to  a  hitherto  undefined  popula¬ 
tion.  The  univariate  program  was  used  to  calculate  measures  of  central 
tendency  and  spread  for  the  Category  X  Cardiovascular  scores.  The 
ANOVA  program  was  used  to  test  for  differences  in  strength,  endurance, 
and  cardiovascular  fitness  among  the  six  test  categories.  In  both 
instances  meaningful  but  guarded  inferences  were  drawn  from  the  test 
data. 

The  capabilities  of  the  TI  59  in  real  world  statistical  analysis 
are  impressive.  The  analyst  can  conduct  sophisticated  analysis  of 
good-sized  samples  unconstrained  by  access  to  large  computers.  Using 
programs  such  as  those  demonstrated  in  this  chapter  the  analyst  need 
not  even  learn  a  programming  language.  All  that  is  required  to  compute 
an  F  ratio  or  prob-value ,  for  example,  is  the  ability  to  follow  simple 
users'  instructions.  While  there  are  certainly  samples  whose  size 
preclude  the  use  of  the  TI  59,  there  are  a  pletheora  of  samples  which 
can  be  analyzed  more  conveniently  and  just  as  efficiently  at  home  or 
at  the  office  using  the  TI  59. 
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TABLE  5-6 


APPENDIX  A 


COMMENTS  ON  SELECTED  LABELS 

LABEL  COMMENT 

A  Directs  action  to  compute  glimpse  probability  of  detection 

for  both  ICWAR  and  I PAR 

|x|  Computes  glimpse  probability  of  detection  for  ICWAR  only 

based  on  target  range 

X  Computes  glimpse  probability  of  detection  for  IPAR  only 

based  on  target  range 

E  and  LnX  Determines  if  target  is  less  than  40  KM  from  unit  and 
directs  processing  accordingly. 

D'  Calculates  target  range  after  elapse  of  target  lock-on  time. 
(SBR)  EXC  and  INT  Increments  range  of  I HAWK  missile  after  firing 
FIX  Prints  24  (no  firing  capability) 

DSZ  Prints  14  (cruise  missile  identified  as  friend) 

NOP  Prints  66  (IFF  is  nonoperational) 

OP  Prints  1  (target  detected)  at  range  (KM) ,  begins  engagement 

sequence 

(SBR)  RAD  and  yx  Increments  simulated  air  battle  time  clock  (R19) . 

i 

(SBR)  GRAD  and  3T  Computes  cruise  missile  rate  of  approach  and  incre¬ 
ments  target  range  (Rqq)  • 

(SBR)  ENG  Random  number  generation  -  normal  distribution 

(SBR)  P— *R  Random  number  generator  -  uniform  distribution 
PRD  Prints  23  (no  detection  capability) 

SIN  Prints  "KILL"  and  range  of  target  kill 

LOG  Prints  "PENETRATION" 

(SBR)  WRITE  AND  TT  Prints  "PENETRATION"  if  target  is  8  KM  or  less 

from  unit. 
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APPENDIX  B 


INSTRUCTION  SHEET 
This  is  a  physical  fitness  test. 

I.  WHAT  WE  ARE  EXAMINING  AND  WHY 

The  following  tests  will  look  at  the  complex  systems  that  make  up 
the  network  of  health  and  fitness.  First  we  will  test  your  upper  body 
endurance  and  muscular  strength.  Next  we  will  examine  your  cardiovas¬ 
cular  health.  This  test  is  a  component  of  a  research  project  being 
conducted  by  graduate  students  of  the  Naval  Postgraduate  School. 

II.  MUSCULAR  ENDURANCE  AND  STRENGTH 

Muscular  endurance  is  often  synonymously  and  incorrectly  used  in 
place  of  muscular  strength.  Muscular  strength  is  the  ability  of  your 
muscular  system  to  exert  maximum  force  against  an  object  or  resistance 
all  at  once,  your  ability  to  exert  a  maximum  force  a  single  time. 

Muscular  endurance  relates  to  the  ability  to  exert  force,  not  necessarily 
maximal,  over  an  extended  time  period.  As  with  all  the  components  of 
fitness,  these  two  concepts  are  interrelated  but  distinctly  different 
from  each  other.  Each  concerns  itself  with  particular  capacities  of 
fitness. 

III.  CARDIOVASCULAR  HEALTH 

Objective:  To  measure  your  heart-rate  response  to  exercise.  This 
test  will  determine  the  relative  efficiency  of  your  heart  and  circula¬ 
tory  system.  Your  heart  rate  can  be  used  to  describe  the  fitness  level 
of  your  body  in  three  areas:  how  much  oxygen  you  need,  how  much  blood 
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your  heart  must  pump  to  supply  this  need,  and  how  hard  your  heart  must 
work  at  this  task.  If  your  need  for  oxygen  is  not  being  fulfilled,  your 
body  is  working  in  an  inefficient  manner.  Consequently,  more  blood  will 
have  to  be  pumped  through  your  circulatory  system  at  a  faster  rate  to 
get  the  oxygen  to  the  muscles  and  organs  that  need  it.  The  heart  has 
the  responsibility  of  satisfying  your  body's  need  for  oxygen.  It  will 
have  to  beat  more  frequently  to  circulate  the  blood  throughout  your 
system.  If  your  body  works  in  an  efficient  manner,  its  need  for 
oxygen  is  being  fulfilled. 

Thus  cardiovascular  health  relates  to  the  ability  of  the  heart, 
lungs,  and  blood  vessels  to  work  in  unison  without  strain.  Regardless 
of  what  the  task  is,  whether  physical  or  mental,  the  cardiovascular 
system  should  be  able  to  handle  it.  When  you  have  high  levels  of 
cardiovascular  health  you  perform  with  more  efficiency  and  you  are 
more  effective  at  what  you  do . 
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QUESTIONNAIRE 


THIS  IS  A  PHYSICAL  FITNESS  TEST 
READ  THIS  DOCUMENT  CAREFULLY  BEFORE  THIS  TEST  3EGINS 

(1)  This  test  is  completely  voluntary!  You  may  decline  testing  now  if 
you  wish.  You  may  stop  at  any  time  during  this  test  and  decline  further 
testing.  You  are  under  no  obligation  to  complete  this  test. 

(2)  Before  proceeding  with  this  test,  you  should  assure  yourself  and 
your  tester  that  there  have  been  no  incidents  in  your  medical  history 
that  would  prohibit  you  from  pursuing  this  testing.  Your  medical  history 
is  relevant  to  this  physical  fitness  test. 


(з)  Please  answer  the  following  statements:  YES  or  NO 

(a)  I  have  a  heart  related  disease. 

(b)  I  have  high  blood  pressure. 

(c)  I  often  feel  faint  and  suffer  spells 
of  dizziness . 

(d)  I  have  recently  or  in  the  past  felt  pain, 
heaviness  or  pressure  in  my  chest. 

(,e)  I  have  felt  pain,  heaviness  or  pressure 
in  my  chest  when  I  walk  uphill . 

(f)  My  doctor  has  advised  me  not  to  engage  in 
physical  exercise  or  physical  activity. 

(и)  I  have  read  and  fully  ’understand  this  document. 

DATE  SIGNATURE 
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DATA  SHEET 


NAME _ 

AGE _  WEIGHT _  DATE _ 

PHYSICAL  ACTIVITY  CATEGORY:  I  II  III  IV  V  VI 

STRENGTH  TEST: 

Maximum  Bench  Press  _ 

Maximum  Curl  _ 

ENDURANCE  TEST: 

Bench  Press  Repetitions  With  100  Pounds 


Curl  Repetitions  With  55  Pounds 


CARDIOVASCULAR  TEST:  (Pulse  Rate) 

Resting 

Immediate  Post  Exercise 
30  Seconds  Post  Exercise 
60  Seconds  Post  Exercise 
120  Seconds  Post  Exercise 
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CATEGORY  I  RESULTS 


SUBJECT 

ACE 

WEIGHT 

CARDIOVASCULAR 
SCORE  CD 

STRENGTH 
SCORE  (2) 

ENDURANCE 
SCORE  (3) 

1 

29 

190 

56.5 

1.26 

24.1 

2 

33 

155 

58.0 

1.32 

12.9 

3 

36 

170 

58.0 

•  97 

5-5 

31 

205 

4 4.5 

1.19 

27.9 

5 

37 

15? 

33-5 

1.13 

10.1 

6 

34 

160 

40.0 

1.16 

11.0 

(1)  See  Table  B-17  for  derivation 

(2)  See  Table  B-18  for  derivation 

(3)  See  Table  B-19  for  derivation 
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CATEGORY  II  RESULTS 


SUBJECT 

AGE 

WEIGHT 

CARDIOVASCULAR 
SCORE  (1) 

STRENGTH 
SCORE  (2) 

1 

33 

16? 

91.5 

1.4? 

2 

34 

195 

96.5 

1M 

3 

35 

180 

55-0 

l.ll 

4 

33 

160 

93-5 

1.44 

5 

28 

127 

82.0 

1.34 

6 

30 

158 

82.0 

1.11 

(1)  See  Table  B-20  for  derivation 

(2)  See  Table  B-21  for  derivation 

(3)  See  Table  B-22  for  derivation 
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ENDURANCE 
SCORE  (3) 


22.0 


28.9 


12.2 


13-6 


6.5 


6.6 


125 


CATEGORY  III  RESULTS 


SUBJECT 

AGE 

WEIGHT 

CARDIOVASCULAR 
SCORE  (1) 

STRENGTH 
SCORE  (2) 

ENDURANCE 
SCORE  (3) 

1 

22 

200 

35.0 

2.46 

50.0 

2 

23 

202 

56.0 

2.35 

62.5 

3 

28 

205 

44.5 

2.49 

37.2 

4 

34 

220 

86.5 

2.23 

42.7 

5 

25 

180 

67.0 

2.83 

54.8 

6 

29 

225 

24.5 

2.00 

jl  .6 

(1)  See  Table  B-23  for  derivation 

(2)  See  Table  B-24  for  derivation 

(3)  See  Table  B-25  for  derivation 
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CATECORY  IV  RESULTS 


SUBJECT 

AGE 

WEIGHT 

CARDIOVASCULAR 
SCORE  (1) 

STRENGTH 
SCORE  (2) 

ENDURANCE 
SCORE  (3) 

1 

28 

200 

96.5 

2.30 

47.5 

2 

31 

180 

79.0 

2.67 

49.0 

3 

32 

165 

85.0 

2.00 

43.7 

4 

32 

160 

79-0 

2.19 

47.2 

5 

19 

240 

89-5 

1.85 

34.3 

6 

26 

180 

71.5 

2.39 

44.6 

(1)  See  Table  B-26  for  derivation 

(2)  See  Table  B-2?  for  derivation 

(3)  See  Table  B-28  for  derivation 
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CATEGORY  V  RESULTS 


SUBJECT 

AGE 

WEIGHT 

CARDIOVASCULAR 
SCORE  (1) 

STRENGTH 
SCORE  (2) 

ENDURANCE 
SCORE  (3) 

1 

30 

165 

39.5 

2.03 

25-5 

2 

24 

175 

55.0 

2.43 

54.1 

3 

19 

185 

47.0 

1.86 

47.0 

u 

34 

179 

77.5 

1.73 

33-5 

5 

35 

205 

59.5 

1.15 

22.0 

6 

30 

120 

56.5 

1.75 

23.8 

(1)  See  Table  B-29  for  derivation 

(2)  See  Table  B-30  for  derivation 

(3)  See  Table  B-31  for  derivation 
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CATEGORY  VI  RESULTS 


SUBJECT 

ACE 

WEIGHT 

CARDIOVASCULAR 
SCORE  (1) 

STRENGTH 
SCORE  (2) 

ENDURANCE 
SCORE  (3) 

1 

27 

180 

61.0 

1.42 

23.1 

2 

32 

125 

^7.5 

1.56 

17.4 

3 

31 

170 

05-3 

1.26 

14.7 

- 

31 

170 

85.0 

1.41 

15-3 

5 

175 

36.5 

1.3^ 

11.4 

6 

30 

205 

30.0 

1.15 

17.9 

(1)  See  Table  3-32  for  derivation 

(2)  See  Table  B-33  for  derivation 
(3/  See  Table  B-34  for  derivation 


TABLE  B-10 


129 


Univariate  User  Instructions 


Step 

Procedure 

Enter 

Press 

Display 

1. 

Initialize 

2nd  E' 

31 

2. 

Enter  data 

X. 

1 

A 

i 

Repeat  for  each  x^ 

3. 

Recall  Statistics 

2nd  2 

X 

INV  2nd  2 

s 

2nd  OP  11 

MSD 

RCL  12 

^in 

RCL  13 

^Snax 

RCL  14 

MIDVAL 

RCL  15 

range 

2nd  A' 

MAD 

RCL  03 

n 

4. 

Recall  data  entered  in 

Step  2  if  desired. 

RCL  31 

xi 

RCL  32 

x2 

RCL  31  +  i 
-1 

xi 
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UNIVARIATE  PROGRAM 
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Enter  data  for  each  row 

A 

ri 

Li 

Reset  pointer  if  more 
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CATEGORY  II  -  ADJUSTED  STRENGTH  SCORES 
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CATEGORY  III  -  ADJUSTED  STRENGTH  SCORES 
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CATEGORY  III  -  ADJUSTED  ENDURANCE  SCORES 
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CATECORY  IV  -  ADJUSTED  STRENGTH  SCORES 


c 

5ENCH  PRESS  + 

CURL 

S4i 

30DY  WEIGHT 

S41 

= 

290  +  170 

200 

2.30 

S42 

305  +  175  _ 

180 

2.67 

S43 

= 

205  +  123 

133 

2.00 

S44 

= 

230  +  120 

160 

2 . 19 

310  ^  I35  _ 

240 

1.85 

S46 

= 

275  +  155 

180 

2.39 

S4x 

270  +  14? 

188 

2.22 

TABLE  3-27 


149 


CATEGORY  IV  -  ADJUSTED  ENDURANCE 


SCORES 


(bodTW)  (NUm*°F  BENCH  PRESS v 
repetitions  )  + 

(30BY~i^  (NnraER  0F  CURD, 

°DY  rfA'  V  repETITIOns  ) 


/lOOv 

200'  (^2)  t  (-J&) 

*-^o  k o0n ; 


2C0^  (76)  =  E7.5 


100,  ,  ,  Cc 

^  k53)  +  (64)  = 


4  \  .  C  r 

1-/v^)+lI (51)  = 


, 100.  ,  , 

<37)  -  ■:-§) 


ToO-)  (?°)  =  47.; 


=  f10Q\  •  ,  rc 

"5  2^0 ;  '"55)  -  (^g)  (50)  = 
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E4x  =  (Jgg)  (49)  +  (-g)  (6l)  =  ^ 
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CATECORY  V  -  ADJUSTED  STRENGTH  SCORES 


S5i 

BENCH  PRESS  + 

CURL 

BODY  WEIGHT 

S51 

= 

220  +  115  _ 
16  3 

2.03 

S52 

= 

265  +  160 

175 

2.L3 

S33 

= 

220  +  125 

185 

1.86 

S54 

= 

200  +  110 

179 

1.73 

S55 

= 

170  +65 

205 

1.15 

S56 

= 

1^0  +70 

120 

1.75 

S5x 

= 

203  +  108 

172 

1.81 
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CATEGORY  V  -  ADJUSTED  ENDURANCE  SCORES 


-  f  100  )  /NUMBER  OF  BENCH  FRESSx  + 

5i  ~  BODY  WT  1  REPETITIONS  1 

/  55  x  , NUMBER  OF  CURLx 

lBODY  WT'  ^  REPETITIONS  1 


E51  = 

/100s 

W 

(30)  ' 

♦<# 

(22) 

=  25-5 

E52  = 

(100) 

475' 

(37) 

(105) 

=  54.: 

E53  = 

(100) 

485' 

(54) 

-  (44 

485' 

(60) 

=  47.0 

E54  = 

(100) 

479; 

(33) 

+  (-55) 
479' 
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=  33-5 

w 

4n 

4r\ 

II 
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V205J 

(28) 

+  ( — a^) 
y205J 
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=  22.0 

E56  = 
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420' 

(22) 

+ 

(12) 

=  23.8 

E5x 

(100) 

4?2; 

(34) 

+  r-55) 
472 ; 

(47) 

=  34.8 
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CATEGORY  VI  -  ADJUSTED  STRENGTH  SCORES 


_  BENCH  PRESS  +  CURL 
b6i  ~  30DY  WEIGHT 


S6l 

= 

155  +  100 
1S0 

S62 

= 

120  +  75 
125 

S63 

= 

130  +  35 
170 

3  64 

= 

135  +  105 
170 

S65 

= 

145  +  90 
175 

S66 

= 

145  +  90 

205 

S6x 

= 

138  +  91 
171 

1.42 

1.56 

1.26 

1.41 

1.34 

1.15 

1.34 
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CATEGORY  VI  -  ADJUSTED  ENDURANCE  SCORES 

E  -  (  100  'l  /NUMEER  OF  BENCH  PRESS, 

oi  ~  ''BODY  WT;  ^  REPETITIONS  ’ 

/  55  x  / NUMBER  OF  CURL, 

lBODY  WT;  ^  REPETITIONS  ‘ 

E6l  =  (25)  +  (30)  =  23,1 

E62  =  ^T§f^  ^  +  ^T253  d4)  =  17,4 

E63  =  (I^}  (14)  +  ^  (2°)  =  14-7 

E64  =  (^§)  (15)  +  (^)  (20)  =  15.3 

E65  =  ("lf§)  (l°)  +  (175)  (i3)  =  11,4 

S66  =  (!§§)  (23)  +  (205)  (25)  =  1^.9 

E6x  =  O  (17)  +  ("171 }  (21)  =  l6’8 
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